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The  study  described  in  this  thesis  was  stimulated 
by  the  need  for  a  working  multipolar  model  of  brain  activity. 
Many  investigators  are  presently  studying  the  magnetic  fields 
emmitted  by  cerebral  cortex  sources,  but  few  of  these  studies 
are  aimed  at  modeling  the  observed  field.  My  strong  interest 
in  biological  systems  of  communication,  computation,  and 
control  has  led  me  to  study  both  biological  system  operation 
and  modern  estimation  and  control  techniques. 

A  degree  of  uncertainty,  caused  by  inaccurate 
measurements,  random  occurances,  and  vague  descriptions  will 
naturally  corrupt  any  attempt  to  estimate  biological  states. 
Fortunately,  these  uncertainties  may  be  reduced  through 
deterministic,  stochastic,  and  fuzzy  variable  methods 
respectively.  A  priliminary  step  required  before  the 
estimation  of  any  real  world  event  is  the  development  of 
proper  models.  Although  these  concepts  have  been  refined  in 
applications  of  flight  control,  interplanetary  navigation, 
etc,  models  of  our  own  cerebral  cortex  are  still  in  their 
infancy.  This  thesis  study  was  therefore  oriented  to  expand 
our  knowledge  in  properly  modeling  the  observed  activity  of 
the  brain. 

The  model  developed  considers  two  or  three  dipoles 
located  in  a  conducting  sphere.  Since  both  a  software  model 
and  an  actual  conducting  sphere  experiment  were  desired, 
time  allowed  only  five  dipole  combinations  to  be 
investigated.  Although  the  study  was  limited  to  only  a  few 


cases,  the  results  are  of  significant  value  and  an  expansion 
of  the  model  should  definitely  be  considered. 
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this  study  helped  me  hurdle  every  threatening  barrier. 
The  list  of  credits  would  certainly  be  incomplete  without 
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environment  and  contributing  to  the  project  without 
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Abstract 

A  Multipole  model  of  magnetic  field  sources  in  the 
cerebral  cortex  was  developed  as  an  expansion  of  Grynszpan's 
single  dipole  in  a  conducting  sphere  model.  Both  computer 
simulation  and  actual  conducting  sphere  experiments  were  used 
to  generate  field  plots  in  five  separate  dipole 
configurations.  These  include  one  case  of  a  single  dipole, 
three  cases  of  three  dipoles  on  a  single  radial,  and  one  case 
of  two  dipoles  on  independent  radials. 

The  results  show  that  dipoles  oriented  on  a  single 
radial  do  not  produce  multiple  flux  peaks,  but  rather  a 
twisted  version  of  the  single  dipole  field.  Sufficiently 
distant  dipoles  on  independent  radials,  however,  do  produce 
multiple  peaks.  The  largest  field  distortions  occurred  when 
nonconducting  masses  were  in  close  proximity  to  the  active 
dipoles  v 
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I.  Introduction 


The  human  central  nervous  system  (CNS)  is  probably  the 
most  complex  operating  network  known  to  man,  capable  of 
manufacturing  its  own  environment  and  even  evaluating  its  own 
sensory  perception.  Besides  higher  intellect,  man  displays 
the  ability  to  consciously  control  precise  muscle  movements, 
a  quality  also  shared  with  other  animals.  In  mammals,  the 
structure  responsible  for  such  movements  is  the  cerebral 
cortex.  Fortunately,  the  cortex  can  be  effectively  studied 
because  of  its  highly  organized  nature  and  the  close 
correspondence  of  its  basic  elements  to  evolutionary 
precursors.  This  characteristic  behavior  of  both  simple  and 
complex  biological  systems  has  contributed  greatly  in 
understanding  at  least  some  aspects  of  the  organic  computing 
process.  Here,  a  review  of  fundamental  neurology  would  prove 
advantageous  in  the  interest  of  developing  a  useful  model  of 
the  observed  electro-magnetic  activity. 


Basic 


(1,2:104-120) 


The  basic  active  element  of  the  nervous  system  is  the 


neuron,  or  nerve  cell.  Unlike  other  cells  in  the  body,  the 
neuron  has  the  ability  to  output  a  sizeable  electrical  pulse 
in  response  to  multiple  inputs.  The  output  pulse,  called  an 
action  potential  (AP) ,  is  the  result  of  an  electro-chemical 
process.  Although  the  neuron  is  a  very  special  cell,  its 


operation  stems  from  basic  electrophysiology  found  in  all 


cells  of  the  body. 


1 


Over  the  millions  of  years  of  evolution,  Earth 
creatures,  including  the  human  animal,  have  retained  a 
fundamental  foundation  of  a  cell  immersed  in  salt  water.  The 
primary  ions  of  interest  are  Sodium  (*3a  +  )  ,  Chloride  (C1-), 
and  Potassium  (K+) .  A  cell's  semi-permeable  membrane  allows 
a  restricted  flow  of  these  ions  from  the  inside  of  the  cell 
to  the  outside,  and  vice  versa.  The  result  is  a  controlled 
potential  difference  across  the  cell's  membrane,  the 
foundation  for  sustained  electrical  activity.  Figure  1-1 
depicts  the  relative  concentrations  of  Na+,  C1-,  and  K+  with 
respect  to  a  neuron. 


Figure  1-1.  Relative  Concentrations  of  Intra  and 
Extra  Cellular  Ions.  The  large  symbols  represent 
a  high  concentration  and  the  arrows  represent  dif¬ 
fusion  . 

Communication  from  neurons  to  other  neurons,  muscle 
fibers,  or  glands  is  primarily  accomplished  via  action 
potentials  (AP) .  With  the  application  of  sufficient 
stimulus,  an  affected  section  of  membrane  experiences  a 
sudden  rise  of  Na+  permeability.  The  rapid  influx  of  Na+ 
results  in  a  localized  breakdown  of  membrane  potential  called 


Figure  II-l.  The  Biot-Savart  Law  (5:233). 


The  magnetic  field  intensity  is  derived  by  considering 
the  generated  force  per  current  element  along  the  path 
length.  The  conducting  length  dL  with  a  current  I  has  a 
"current  moment"  given  by:  IdL  =  current  x  length  =  current 
moment.  The  Magnetic  Flux  Density  B  is  defined  as  the  force 
per  current  moment  (14:146):  B  =  dF/ldL.  The  magnetic  flux 
density  thus  has  the  units  of  newton  per  ampere-meter.  These 
units  are  equivalent  to  the  units  weber  per  meter  squared, 
10,000  gauss,  or  the  International  System  of  units,  Tesla 
(5:260) . 

The  magnetic  flux  density  conceptually  refers  to  the 
number  of  magnetic  flux  lines  passing  through  a  unit  area. 
It  is  related  to  the  magnetic  field  intensity  through  the 
properties  of  the  medium  involved  (7:29): 
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II.  Field  Measurement  Theory 

Ms±£  Electromagnetics 

The  theory  relating  to  magnetic  fields  generated  by 
current  flow  is  based  upon  the  Biot-Savart  law.  As  noted  in 
the  introduction,  this  law  states  that  at  any  measuring 
point,  the  magnitude  of  the  magnetic  field  intensity  is 
proportional  to  the  dipole  current,  length,  and  the  angle 
between  the  dipole  and  the  measuring  point  (5:232).  To  avoid 
any  generalities,  the  exact  law  should  include  differential 
dipole  lengths  and  vector  parameters.  Equation  (1)  describes 
the  magnetic  field  intensity  generated  by  a  differential 
dipole  length: 


IdL  x  a 

dH  =  -  {1) 

47rR2 


where 

dH  =  infinitesimal  magnetic  field  intensity 
dL  =  infinitesimal  path  length 
I  =  current 

a  =  unit  vector  from  path  to  measuring  point 
R  =  distance  between  path  and  measuring  point 

("x"  denotes  the  vector  cross  product) 


Figure  II-l  depicts  how  the  differential  length  of  equation 
(1)  contributes  to  the  magnetic  field  intensity  measured  at 
point  P.  Integration  of  dH  along  the  entire  current  path 
yields  the  generated  field. 


15 


The  remaining  five  chapters  of  this  report  will  cover 
existing  magnetic  field  measurement  theory,  the  multipole 
model  developement ,  experimental  procedure,  discussion,  and 
conclusions . 

Chapter  II  includes  some  basic  electromagnetics  and 
fundamentals  in  SQUID  technology. 

In  chapter  III  a  multipole  model  is  developed  as  a 
direct  extension  of  Grynszpan's  single  dipole  model.  The 
method  of  using  meridians  and  parallels  to  navigate  the 
sphere  is  also  explained,  as  well  as  the  software  used  to 
develop  magnetic  contour  maps. 

Chapter  IV  is  a  detailed  discription  of  the  experimental 
equipment  and  procedure.  This  includes  the  particular  2nd 
order  SQUID  gradiometer  used,  the  gimbaled  sphere  mount, 
support  equipment,  and  the  procedure  of  collecting  data. 

The  results  are  presented  in  chapter  V.  A  comparison  is 
made  between  the  theoretical  and  the  experimental  plots. 

The  discussion  is  given  in  chapter  VI.  A  suggestion  of 
how  an  array  of  such  dipole  stacks  can  be  used  as  an 
estimator  is  also  presented.  Finally,  the  conclusion  will 
include  an  overview  of  the  project  and  some  recommendations 
for  continuing  studies. 


Scope  And 


The  multipole  model  developed  consists  of  three  current 
dipoles  within  a  saline  filled  pyrex  sphere.  Three  primary 
assumptions  exist  in  this  model.  First,  the  skull  is  assumed 
to  closely  resemble  a  non-conducting  spherical  shell. 
Secondly,  this  spherical  shell  contains  a  homogeneous 
conducting  fluid.  And  thirdly,  the  source  of  interest 
consists  of  a  tight  group  of  localized  dipoles.  The  first 
two  assumptions  eliminate  the  complexity  associated  with  the 
true  shape  and  consistency  of  the  head.  Although  the  head  is 
not  perfectly  spherical,  many  investigators  accept  the 
spherical  model  as  reasonable  for  most  areas  of  interest 
(10:27).  The  inhomogeneties  of  the  midbrain,  cerebral 
cortex,  cerebral  spinal  fluid,  meninges,  skull,  etc,  will 
distort  and  attenuate  the  detected  field.  Reputable 
investigators,  however,  claim  that  these  effects  do  not 
seriously  disrupt  source  localization  (10:27). 

Once  a  source  is  localized,  a  multipole  model  can  be 
adjusted  to  best  represent  the  distortions  found  in  the 
observed  field.  Two  cases  are  investigated.  In  the  first, 
the  dipoles  comprise  a  single  stack  alligned  on  a  radial  and 
are  oriented  in  variable  directions  perpendicular  to  the 
radial.  The  second  case  considers  tangential  dipoles  lying 
on  different  radials.  A  comparison  will  be  made  between 
computer  generated  isofield  maps  and  actual  isofield  maps 
detected  by  a  second  order  SQUID  gradiometer. 


application  of  the  effecting  stimulus.  The  model  developed 
in  this  thesis,  however,  represents  an  isolated  source  of 
cortical  activity.  Whereas  averaging  is  currently  necessary 
in  actual  brain  studies,  it  will  not  be  necessary  in  a  model 
with  deterministic  inputs. 

The  problem  statement  motivating  the  model  will  be 
discussed  next  and  a  detailed  description  of  the  theory  and 
model  developement  will  follow  in  subsequent  chapters. 

Bijpblfim  Statement 

Many  researchers  of  today  who  are  studying  the 
characteristics  of  evoked  responses  in  the  cerebral  cortex 
isolate  the  source  of  activity  by  assuming  a  model  of  a 
single  current  dipole.  The  dipole  is  located  between  the 
positive  and  negative  flux  density  peaks  and  at  an  easily 
determined  depth  and  orientation  (7:130-131).  However,  the 
observed  two  dimensional  isofield  maps  resemble  distorted 
dipole  fields  and  usually  exhibit  several  flux  density  peaks. 
The  departure  from  circular  fields  indicates  various  boundary 
conditions  and  together  with  multiple  peaks  suggests  that 
more  than  just  one  dipole  should  model  the  isolated  source. 
In  1983,  Katila  and  Karp  published  a  multipole  expansion  of 
cardiac  magnetic  fields  seeking  to  find  the  "simplest 
description  of  an  equivalent  source  that  accounts  for  the 
observed  field  pattern"  (7:252-254).  This  thesis  will  apply 
a  similar  thought  process  in  modeling  the  observed  field  of 
isolated  cerebral  cortex  sources. 
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locked  responses,  a  technique  also  used  on  EEG  and  MEG  data 
to  eliminate  noise  and  localize  dipole  sources. 

Each  of  the  three  data  aquisition  methods  mentioned  have 
advantages  and  disadvantages,  and  in  the  interest  of  studying 
the  human  cerebral  cortex,  these  qualities  must  be  evaluated. 
The  microelectrode  histograms  can  be  eliminated  immediately 
since  real  time  analysis  is  impossible  and  sticking 
electrodes  into  the  cortex  of  human  subjects  is  out  of  the 
question.  Although  the  pros  and  cons  between  magneto¬ 
encephalography  and  electro-encephalography  deserve  more 
attention,  the  particular  application  of  this  thesis  favor 
the  MEG's  advantages  over  those  of  the  EEG.  Basicly,  the 
four  deciding  characteristics  are: 

(1)  The  SQUID  gradiometer  detects  only  field  compo¬ 
nents  of  sources  close  to  the  sensor  (9:230). 

(2)  The  predominant  fields  sensed  are  due  to  the 
intracellular  currents  with  little  effect  from 
the  extracellular  volume  currents  (9:237). 

(3)  The  magnetic  field  is  not  affected  by  insulating 
layers,  which  minimizes  the  effects  of  the  non¬ 
computing  meninges,  skull,  and  skin. 

(4)  The  MEG  is  a  non-invasive  method  of  obtaining  brain 
waves  without  the  added  complexity  of  reference 
electrodes . 

From  this  point  forward,  reference  to  data  aquisition  will  be 
in  terms  of  magneto-encephalography. 

The  best  known  method  of  determining  the  location  of  a 
dipole  source  and  observing  its  activity  is  through 
averaging.  Because  of  the  physical  characteristics  of 
chemical  transmitters  and  propagating  action  potentials,  at 
least  some  neural  responses  are  time-locked  to  the 


The  recording  of  continuous  magnetic  fields  above  the 
surface  of  the  head  is  termed  magneto-encephalography  (MEG) . 
Since  the  first  known  MEG  recording  by  David  Cohen  in  1967  at 
the  Massachusetts  Institute  of  Technology  (6:784-786),  the 
magneto-encephalography  technique  has  been  substantially 
improved.  The  present  day  MEG  is  obtained  using  a 
Superconducting  Quantum  Interference  Device,  commonly 
referred  to  as  a  SQUID.  In  the  interest  of  eliminating 
distant  noises  of  large  magnitude,  contemporary  SQUIDs  are 
constructed  as  first  and  second  order  gradiometers.  The 
multiloop  gradiometer  operating  at  approximately  4.2  degrees 
Kelvin  (7:69)  allows  the  SQUID  to  detect  extremely  weak 
magnetic  fields  (on  the  order  of  10exp(-15)  Tesla) 
uncorrupted  by  large  fields  such  as  room  lights,  power 
supplies,  etc.  The  penalty  paid  for  the  reduction  in  noise 
is  a  reduction  in  sensitivity  for  each  additional  flux 
detecting  loop.  This  penalty  is  considered  relatively  small 
when  there  is  a  need  to  operarate  a  SQUID  in  the  noisy 
environments  of  unshielded  laboratories  and  hospitals,  hence 
the  use  of  a  "second  order  gradiometer"  is  advised  (7:92-93). 
Concepts  about  SQUID  and  second  order  gradiometer  operation 
will  be  covered  in  chapters  II  and  III  respectively. 

As  discussed  above,  the  low  frequency  waves  containing 
information  about  the  on-going  brain  activity  can  be  detected 
through  EEG  or  MEG  methods.  Frank  R.  Ervin  has  even  detected 
these  same  low  frequency  "brain  waves"  by  constructing 
histograms  of  microelectrode  data  from  single  cortical 
neurons  (8:35-51).  Ervin's  data  is  the  average  of  many  time- 


10 


electric  fields"  (4:21).  The  resultant  summation  wave  of 
this  electrical  activity  yields  a  view  of  the  general  brain 
state.  Interpretation  of  the  EEG  requires  a  thorough 
evaluation  of  the  wave  frequency,  amplitude,  form, 
periodicity,  and  also  of  the  electrode  locations.  The 
detected  brain  wave  frequencies  are  relatively  low  and  are 
classified  into  four  broad  bands  (4:50): 


alfa - 8  through  13  hertz 

beta - 14  and  above  hertz 

theta - 4  through  8  hertz 

delta - 0.5  through  3  hertz 


The  amplitude  of  scalp  potentials  is  on  the  order  of 
microvolts  and  evaluated  as  relative  peaks  and  troughs. 
Patterns  described  by  the  above  criteria  have  been  studied  in 
both  normal  and  abnormal  brains  and  form  the  basis  for 
clinical  determination  of  neural  abnormalities.  Reference 
4:131-236  for  a  clinical  evaluation  of  various  EEG  patterns. 

Electrical  brain  activity  and  current  flow  also  create 
distinctive  magnetic  fields  reflecting  the  source  of 
activity.  If  the  sources  are  modelled  as  simple  dipoles  (two 
equal  but  oppositely  charged  particles  separated  in  space) , 
then  the  detected  magnetic  field  is  directly  related  to  the 
dipole  source  current  by  the  Biot-Savart  law.  This  law, 
first  proposed  in  1820,  will  be  discussed  further  in 
subsequent  chapters.  In  short,  the  law  states  that  at  any 
measuring  point,  the  magnitude  of  the  magnetic  field 
intensity  is  proportional  to  the  dipole  current,  length,  and 
the  angle  between  an  equivalent  dipole  and  the  measuring 


dipoles  contribute  significantly  to  the  magnetic  field 
outside  the  head,  assuming  the  head  can  be  approximated  as  a 
conducting  sphere.  This  result  stems  from  Grynszpans  dipole 
model  which  will  be  discussed  in  chapter  III.  Note  that 
magnetic  field  studies  will  be  expected  to  localize  cerebral 
activity  in  the  columns  found  in  the  cortical  folds. 

Although  the  cerebral  cortex's  general  use  of  the 
columns  is  still  being  debated  in  scientific  circles,  the 
most  promising  theory  is  that  of  two  dimensional  Fourier 
analysis  described  by  Kabrisky  in  1966  (18) .  However,  before 
the  actual  mathematical  computations  performed  by  the  network 
of  cortical  columns  can  be  studied  non-evasively  through 
epidermal  detectors,  the  instrumentation  and  models  for  such 
detectors  must  be  refined. 

Encephalfigr. aphy 

Electrical  activity  implies  current  flow.  Recall  that 
the  CNS  has  been  shown  to  have  a  communications  network  based 
upon  a  transmembrane  flow  of  ions  and  propagating  action 
potentials.  With  the  intense  activity  of  propagating  action 
potentials,  an  induced  current  flows  throughout  the  body 
fluid.  The  induced  current  can  be  subdivided  into  the  two 
classes  of  (1)  intracellular  currents  and  (2)  extracellular 
volume  currents.  The  well  known  electro-encephalogram  (EEG) 
provides  information  about  brain  activity  by  detecting 


surface  potentials  created  from  induced  extracellular  volume 
currents.  In  1853  H.  Helmholtz  noted  that  a  detected  surface 
potential  may  represent  an  "infinite  variety  of  internal 


The  axon  of  most  pyramidal  cells  extends  out  of  the  cortex  to 


connect  to  other  cortex  areas  or  subcortical  areas.  Thus, 
while  the  stellate  cells  are  primarily  cortical 
interconnecting  neurons,  the  pyramidal  cells  act  as  the 
primary  output  neurons. 

Anatomists  have  defined  six  cortical  layers  by  the 
density  of  various  occuppying  cell  types.  However,  some 
areas  such  as  the  motor  cortex  have  many  vertical  pyramidal 
cells  and  few  stellate  or  other  type  cells.  Similarily, 
other  areas  such  as  the  sensory  cortex,  are  primarily 
horizontal  interconnecting  stellate  cells.  Although  the  six 
layers  in  these  cortical  regions  are  hard  to  distinguish, 
their  near  homogeneous  anatomical  structures  are  an  aid  to 
electromagnetic  source  modeling.  Figure  1-5  shows  a  section 
of  mouse  cortex  stained  by  the  Golgi  technique.  Note  the 
highly  vertical  orientation  in  this  region  of  cortex. 


Figure  1-5.  Cerebral  Cortex.  (From  Lorente  de 
No';  Reprinted  in  reference  16:263). 
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Figure  1-4.  A  Neuron  witn  Many  Synaptic 
Junctions  (3:247) 


lira  Cerebral  Cox  tax 

The  cerebral  cortex  is  the  most  recent  addition  to  the 
complex  neurological  systems  found  in  mammals.  In  humans, 
this  sheet  of  neurons,  interconnections,  and  glial  cells 
is  about  2.5  square  feet  in  area  and  about  2  or  3  millimeters 
thick. 

Because  of  its  late  apearance  in  the  evolutionary  chain, 
the  primary  areas  of  the  cortex  are  called  the  neocortical 
cortex.  The  neocortex  consist  of  two  principle  cell  types: 


1.  Stellate  cells. 

2.  Pyramidal  cells, 


The  stellate  cells  are  small  and  have  short  axons  that  stay 
within  the  local  cortex  area.  The  pyramidal  cells,  however, 
have  long  "apical"  dendrites  that  extend  to  the  surface  of 
the  cortex,  and  "basal"  dendrites  that  spread  horizontally. 
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Figure  1-3.  Synaptic  Junctions  between  Neurons 


ix 


depolarization.  Within  milliseconds  the  membrane  is 
repolarized  by  a  lowering  of  Na+  permeability  and  an  increase 
of  K+  permeability.  The  activity  at  one  section  of  membrane 
provides  sufficient  stimulus  to  trigger  adjacent  membrane 
sections.  Thus,  the  AP  propagates  the  length  of  the  neuron 
without  attenuation. 

Altnough  many  inputs  may  be  applied,  each  neuron  is 
equipped  with  only  one  output  channel,  called  an  "axon".  The 
terminating  end  of  the  axon  is  called  a  "synapse"  and  forms  a 
synaptic  junction  with  the  receiving  cell.  An  axon  may  have 
more  than  one  synapse.  Upon  receipt  of  an  AP,  vesicles 
within  the  synapse  release  chemical  transmitters  which  excite 
(produce  an  EPSP-excitatory  post  synaptic  potential)  or 
inhibit  (produce  an  IPSP-inhibitory  post  synaptic  potential) 
the  receiving  cell.  A  multitude  of  EPSPs  and  IPSPs  will 
summate  both  temporally  and  spatially  within  a  receiving 
neuron.  The  neuron  will  then  fire  (initiate  an  AP)  only  if 
the  graded  summation  exceeds  a  fixed  threshold.  Thus,  a 
system  for  computation  and  communication  can  be  achieved  via 
a  network  of  functional  neurons.  A  propagating  AP  is 
diagrammed  in  figure  1-2.  Figure  1-3  depicts  presynaptic 
neurons  synapsing  with  postsynaptic  neurons.  Figure  1-4 
shows  how  a  neuron  may  have  a  multitude  of  inputs,  but  only 


one  output. 
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B  =  u  (H  +  M) 

O 

where 

H  =  magnetic  field  intensity 
L  =  path  length 

M  =  magr etization :  the  response  of  a  material 
wnen  subject  to  a  magnetic  field. 

u.3  =  permeability  of  free  space  (47 r  xlO  7  Hen/m) 

Note  from  equation  (2)  that  the  mere  presence  of  a  material 
causes  an  increase  in  flux  density  over  that  of  free  space. 

This  thesis  is  only  concerned  with  the  single  radial 
componemt  of  the  magnetic  field  vector.  In  this  case,  where 
B  and  M  are  parallel  to  H,  equation  (2)  reduces  to  its  scaler 
form.  A  further  simplif ication  can  be  made  if  the  medium 
involved  is  considered  to  have  properties  similar  to  that  of 
free  space,  i.e,  M=0. 


I  squid  Teghnology 

When  the  interest  of  measuring  extremely  small  magnetic 
fields  developed,  investigators  discovered  that  the  quantum 
effects  of  superconducting  currents  is  related  to  the 
magnetic  flux  density  passing  through  a  superconducting 
current  loop.  This  phenomenon  forms  the  basis  of  the  SQUID 
gradiometer  described  in  chapter  I.  An  overview  of  the 
concepts  will  be  given  in  the  following  paragraphs.  A  more 
indepth  discussion  can  be  found  in  chapter  4  of  reference  7. 

At  temperatures  close  to  absolute  zero,  as  required  for 
superconductivity,  the  attractive  force  between  electrons 
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becomes  more  important  than  disordering  thermal  effects.  As 
a  result,  the  electrons  bind  together  in  pairs  which  appear 
to  have  the  same  quantum  state.  With  a  large  number  of 
electrons  behaving  in  the  same  manner,  a  superconducting 
metal  exhibits  some  unusual  properties.  The  quantum  effect 
on  the  wavelike  properties  of  the  superconducting  electrons 
made  the  development  of  the  SQUID  possible. 

In  1960  Giaever  discovered  a  tunneling  effect  of 
superconducting  electrons  through  an  insulating  barrier 
separating  a  superconductor  and  a  normal  metal  or  two 
superconductors  (15:3,4).  Because  of  the  quantum  behavior  of 
the  electrons,  their  wavelike  properties  can  not  drop 
abruptly  to  zero  at  the  barrier.  Such  an  abrupt  drop  of  wave 
amplitude,  according  to  quantum  theory,  would  cause  a 
gigantic  particle  kinetic  energy.  Instead,  the  "electron 
wave-like  properties"  penetrate  the  barrier  and  decay 
smoothly  to  zero.  Penetration  of  the  barrier  from  both  sides 
allows  the  electrons  to  "tunnel"  from  one  side  to  the  other. 
In  1962  Josephson  predicted  an  interesting  occurance  if  the 
two  superconductors  were  separated  by  a  very  thin  insulating 
barrier.  The  tunneling  electrons,  according  to  Josephson, 
would  set  up  a  supercurrent  through  the  insulator,  i.e, 
current  would  flow  through  the  insulator  without  any  voltage 
drop.  A  superconducting  barrier  is  thus  called  a  Josephson 
j  unction . 
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Loop  with  a 
weak  link 


Equivalent 

circuit 


Figure  II-2.  A  Superconducting  Loop  with 
a  Josephson  Junction  Weak  Link. 


In  addition,  Josephson  showed  that  the  superconducting 


current  is  described  by  (7:74): 


where 


I  =  I  sin 

c 


=  superconducting  current 

=  critical  current  value:  If  the  junction  is 
forced  to  carry  a  current  greater  than  I  , 
a  voltage  drop  must  appear.  This  additional 
current  flows  across  the  junction  in  the 
normal  fashion  (with  resistance) . 

=  phase  difference  between  the  electron  wave 
functions  on  either  side  of  the  junction. 


At  the  Bell  Telephone  Laboratories  Rowell  and  Anderson  were 
the  first  to  experimentally  varify  this  theory  (15:66).  The 
phase  difference  <t>  is  free  to  adjust  itself  to  allow  any 
current  to  flow  up  to  the  critical  current  value.  This 

property  is  called  the  "DC  Josephson  effect"  and  equation  (3) 
is  called  the  "DC  Josephson  equation". 
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If  a  voltage  V  is  present  across  the  Josephson  junction, 
it  turns  out  to  be  proportional  to  the  rate  of  change  of  the 
electron  waves'  phase  difference  0  (7:74): 


where 


h 

V  = - 

4  e 


dt 


2  dt 


(4) 


h  =  Plank's  constant  =  6.63  x  10~34  J*s 
2e  =  charge  on  the  electron  pair 
<J>0  =  flux  quantum  =  h/2e  =  2.07  x  10-7  Weber 


Note  that  a  voltage  appearing  across  the  junction  will 
cause  the  phase  difference  to  change  with  time.  Putting  this 
into  equation  (2)  sets  up  an  alternating  current  and  is  thus 
called  the  "AC  Josephson  effect". 

When  a  flux  density  is  applied  to  a  superconducting  loop 
containing  a  Josephson  juction,  the  induced  voltage  alters 
the  properties  of  the  junction  and  a  non-linear  relationship 
is  established  between  the  applied  magnetic  field  and  the 
current  flowing  around  the  loop.  The  sensitivity  of  a 
Josephson  junction  will  be  on  the  order  of  a  flux  quantum. 

There  are  two  primary  methods  of  measuring  the  effects 
the  applied  flux  has  on  the  Josephson  junction  loop.  Both 
methods  use  the  instrumentation  concept  of  "Parametric 
Amplification".  The  signal  detecting  circuit  is  driven  by  a 
"pump"  voltage  or  current  source.  How  the  pump  source 


affects  the  detecting  circuit  and  vice-versa  is  governed  by  a 
reactive  component.  The  method  of  signal  detection  is  based 
on  the  fact  that  the  capacitance  or  inductance  of  this 
component  is  controlled  by  the  input  signal. 

The  first  of  the  two  systems  is  the  "DC-SQUID".  This 
method  requires  only  a  DC  pump  current  since  the  reactance 
affecting  AC  current  is  produced  by  the  "AC  Josephson 
effect".  A  circuit  diagram  of  DC-SQUID  is  shown  in  figure 
II-3.  Note  that  two  Josephson  junctions  are  required.  If 
just  one  junction  was  used,  the  current  (which  must  be  larger 
than  the  critical  current  to  assure  the  AC  Josephson  effect) 
would  take  the  path  of  least  resistance  and  flow  around  the 
junction. 


Figure  II-3.  Circuit  Diagram  of  a  DC-SQUID 


If  the  pump  current  is  greater  than  the  sum  of  the  two 
critical  currents,  then  both  junctions  will  be  driven  into 
the  "AC  Josephson"  state.  An  input  flux  to  the  coil  will 
alter  the  SQUID  loop  inductance  which  in  turn  alters  the 
voltage  V  measured  across  the  SQUID. 

The  more  commonly  used  flux  detection  system  is  the  RF- 
SQUID .  Unlike  the  the  DC  system  described  above,  the  RF- 
SQUID  requires  only  one  Josephson  junction.  A  circuit 
diagram  is  shown  in  figure  II-4. 


Figure  II-4.  Circuit  Diagram  of  a  RF-SQUID 
(7:77)  . 


In  the  RF  method,  the  loop  inductance  of  the  SQUID 
circuit  is  inductively  coupled  to  a  resonant  tank  circuit 
which  is  driven  by  a  radio  frequency  pumping  current. 
After  resonation  is  set  up  in  the  resonant  tank  circuit,  the 
current  flowing  between  the  capacitor  and  inductor  (see 
figure  II-4)  will  exceed  the  RF  pumping  current  by  the 


quality  factor  "Q".  The  resulting  voltage  across  the 
parallel  capacitor-inductor  circuit  is  amplified  and 
rectified  to  give  the  measured  output  voltage. 

Since  the  SQUID  loop  will  oscillate  at  the  "Josephson 
frequency"  given  by  f  =  2e/h  =  483  MHz/uV  (7:74),  inductive 
coupling  is  best  achieved  at  radio  frequencies.  The 
coupling  between  the  tank  and  SQUID  circuits  allows  the  RF 
pumping  current  to  drive  the  Josephson  junction  into  its  AC 
Josephson  state.  When  flux  density  is  applied  to  the  SQUID 
loop,  the  SQUID  inductance  is  altered,  which  then  alters  the 
resonant  circuit  coupled  to  it.  Changes  in  the  resonant  tank 
circuit  are  then  detected  at  the  output  voltage  V  . 

In  both  the  DC  and  the  RF  systems,  the  output  voltage  is 
a  non-linear  function  of  input  flux  density,  and  periodic 
with  the  period  given  by  the  flux  quantum.  Since  it  is 
desireable  to  compensate  for  the  non-linearities  and  avoid 
the  periodicity  of  the  SQUID  output,  investigators  have 
designed  a  simple  closed  feedback  loop,  from  which  the 
nulling  signal  is  picked  off  as  the  output  measurement. 
Because  the  second  order  gradiometer  used  in  this  thesis 
consists  of  a  RF-SQUID  and  the  feedback  loop  described  above, 
the  combined  system  will  be  described  further. 

Recall  the  RF-SQUID  shown  in  figure  II-4.  Figure  II-5 


shows  the  same  circuit,  only  this  time  it  is  embedded  in  a 
flux  feedback  loop. 


Figure  II-5.  SQUID  Feedback  Loop  (7:81). 


To  gain  an  understanding  of  the  operation  of  the 
feedback  loop,  one  can  start  at  the  feedback  resistor  R  and 
go  around  the  loop.  On  the  feedback  line  there  will  be  a 
modulation  frequency  that  originated  from  the  audio  frequency 
generator  AF  .  This  modulation  frequency  introduces  a  flux 
to  the  SQUID  loop.  Since  the  SQUID  signal  will  be  riding  on 
the  AF,  the  AF  must  be  chosen  to  have  an  appropiately  large 
enough  bandwidth,  usually  about  50KHz  (7:81).  After  passing 
through  an  amplifier  and  a  RF-amplitude  modulation 
detector (diode) ,  the  audio  frequency  modulated  signal  goes 
into  a  "lock  in"  amplifier  which  is  locked  in  on  the  audio 
frequency.  The  "lock  in"  amplifier  simply  consists  of  a 
Phase-Sensitive  Detector  and  a  low  pass  filter.  The 
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modulated  signal,  which  represents  the  SQUID  detected  input 


flux,  then  passes  through  the  feedback  resistor.  The  coil 

after  the  feedback  resistor  is  wound  such  that  the  flux 

signal  introduced  cancels  with  the  flux  present  in  the  SQUID 

loop.  The  negative  feedback  loop  is  thus  closed.  The  null  - 

current  can  be  detected  by  measuring  the  voltage  across  the 

feedback  resistor.  When  there  is  no  input  flux  to  the  SQUID 

coil,  the  negative  feedback  signal  perfectly  cancels  the 

present  SQUID  flux,  and  the  null  current  is  zero.  Similarly, 

when  an  input  flux  is  applied,  the  feedback  system  attempts 

to  compensate  for  it,  and  a  null  current  is  generated.  It 

turns  out  that  the  null  current  is  linearly  related  to  the 

input  flux.  The  system  is  called  a  "flux  locked  loop". 

Flux  EolaLity  and  Sinusoidal  Eh&££  Shift  ‘  ’ 

Most  investigative  studies  use  a  sinusoidal  input  to 
stimulate  cerebral  or  mid  brain  dipoles  of  interest.  At  any 
instant  of  time,  magnetic  flux  lines  will  exit  the  scalp 
surface  on  one  side  of  the  dipole  and  enter  on  the  other. 

See  figure  II-6.  This  is  in  accordance  with  the  right  hand  ( 

rule,  which  lets  the  right  hand  thumb  represent  current 
direction  and  the  curl  of  the  fingers  represent  magnetic  flux 
flow.  The  accepted  convention  lets  the  exiting  (upward)  flux 
be  positive  and  the  entering  (downward)  flux  be  negative. 

Assuming  no  transmission  line  delays  or  reactance,  the 
generated  dipole  flux  will  track  the  input  voltage. 

* 
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P  =  Dipole  moment,  B  =  Magnetic  flux  density, 
J  =  induced  extracellular  volume  current. 


Recall  that  the  SQUID  detection  system  depends  on  flux 
linkages.  The  electromotive  force  (emf)  generated  by 
magnetic  flux  passing  through  a  coil  is  given  by  Faraday's 
law  of  induction: 

6<P 

emf  =  -  N  -  (5) 

dt 


where , 


0  =  Magnetic  flux 
N  =  #  of  turns  in  the  coil. 


Note  that  a  sinusoidal  input  flux  to  a  coil  will  result  in  a 
sinusoidal  emf  with  a  90  degree  phase  shift.  Using 
equation  (5)  and  assuming  negligible  flux  leakage,  it  can  be 
shown  that  the  voltage  (emf)  across  a  secondary  transformer 
coil  is  directly  proportional  to  the  voltage  across  the 
primary  coil: 
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where , 


V2 


v2 =  voltage  across  secondary  coil 
vL  =  voltage  across  primary  coil 
N, =  #  of  turns  in  first  coil 
N, =  #  of  turns  in  second  coil 


Figure  II-7  illustrates  the  inductive  coupling  and 
phase  relations  found  in  a  RF  SQUID. 


Figure  II-7.  Phase  Relation  Between  Dipole  Voltage  and 
Induced  SQUID  Electromotive  Force  (emf). 
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A  sinusoidal  flux  input  to  the  detection  coil  results  in 
a  90  degree  lagging  sinusoid  induced  across  the  input  coil 
and  cor respondently  in  the  SQUID  loop.  The  signal  is 
embedded  in  the  nonlinear  characteristics  of  the  SQUID,  and 
then  relinearized  through  the  feedback  SQUID  control  circuit. 
Although  the  induced  emf  in  the  SQUID  loop  lags  the  input 
flux  by  a  known  90  degrees,  reactances  within  the  rest  of  the 
control  circuitry,  signal  processing  amplifiers,  and  filters 
will  cause  additional  phase  shifts. 

Since  the  reactances  will  shift  sinusoids  of  a  given 
frequency  by  the  same  amount,  the  180  degree  phase  difference 
between  the  +Flux  induced  emf  and  the  -Flux  induced  emf  is 
preserved.  This  character istic  is  used  to  determine  the 
detected  flux's  polarity  by  recording  the  phase  difference 
between  the  input  stimulus  and  the  processed  output  signal. 
The  output  associated  with  +Flux  will  be  180  degrees  out  of 
phase  with  the  output  associated  with  -Flux. 
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III.  Model  Development  and  Plotting 


Expansion  ol  Grynszpan1 s  Model 

In  a  dissertation  published  in  1971  (11:26-32), 
Grynszpan  developed  equations  describing  the  magnetic  field 
produced  by  a  single  dipole  in  a  sphere.  Cuffin  and  Cohen 
outlined  these  results  while  evaluating  the  magnetic  fields 
outside  various  volume  conductor  shapes  (12:374-375).  The 
equation  of  interest  is  the  radial  component  of  the  magnetic 
field  strength.  In  spherical  coordinants,  this  component  is 
the  magnitude  of  the  magnetic  field  vector,  and  hence  the 
magnetic  strength  detected  by  a  SQUID. 

Grynszpan  used  a  spherical  coordinant  system  with  the 
origin  coinciding  with  the  sphere's  origin.  The  dipole  was 
fixed  at  X=Y=0  and  variable  on  the  z-axis.  The  dipole 
orientation  was  fixed  in  the  positive  X  direction  and 
parallel  to  the  XY-plane.  Figure  III-l  depicts  this 
arrangement.  The  calculated  radial  component  is  given  by: 


P  a  Sin  0  Sin  0 


H  = 


47T  R2  y 


2  V 


where 


2  a  Cos  0 


Dipole  Moment 

Dipole  location  on  the  z-axis 
3.14159 

Field  Location  in  Spherical  Coordinants 
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Figure  III-l.  Grynszpan's  Dipole  Model.  The  dipole 
is  represented  by  the  heavy  arrow  at  z=a  and  H(9,^f,R) 
is  the  calculated  radial  field  component. 


Note  that  Grynszpan's  model  considers  only  one 
tangential  dipole  which  has  strict  orientation  constraints. 
It  also  includes  all  derivatives  in  the  Taylor  series 
expansion.  Since  the  multipole  model  intended  for  this 
thesis  requires  various  tangential  dipole  positions,  and  the 
2nd  order  SQUID  gradiometer  does  not  detect  the  Oth  and  1st 
order  derivatives,  certain  modifications  were  made  to  the 
calculated  radial  field  component. 


A  super  insulated  fiberglass  dewar  houses  the  BMP  probe 
and  its  liquid  helium  bath.  The  gradiometer  coils  are 
positioned  in  a  thin  neck  pertruding  from  the  bottom  of  the 
dewar.  This  allows  the  gradiometer  to  be  placed  close  to  the 
magnetic  field  source. 

Recall  the  RF  SQUID  and  feedback  system  described  in 
chapter  II  (see  figure  1 1  —  5 ) .  The  S.H.E  system  is  a  similar 
RF  SQUID  with  the  exception  that  the  feedback  path  is  fed 
directly  into  the  BMP  probe  output  instead  of  being 
inductively  coupled  to  the  SQUID  loop.  Conceptual  circuit 
operation,  however,  is  the  same.  The  output  of  the  BMP  probe 
is  fed  into  a  S.H.E.  model  300X  RF  head  unit  mounted  on  the 
top  of  the  dewar.  The  RF  head  connects  to  the  S.H.E.  model 
30X  control  unit.  Together  the  system  is  called  the  S.H.E 
330X  and  is  diagrammed  in  figure  IV-2. 


Figure  IV-1.  S.H.E.  Model  BMP  Biomagnetic  Probe 
(Reference  19  and  20). 


IV.  The  Multipole  in  a  Sphere  Experiment 


Five  dipole  configurations  were  examined.  The  first  run 
consisted  of  a  single  dipole  located  on  the  z-axis  and 
oriented  in  the  positive  x  direction  as  described  by 
Grynszpan's  original  model  (see  figure  III-l).  This 
configuration  povided  a  reference  for  the  behavior  of  a 
single  dipole  within  a  conducting  sphere.  The  second  run 
extended  to  a  three  dipole  stack  on  the  z-axis.  The  dipoles 
were  all  oriented  in  the  positive  x  direction  but  located  at 
different  z-axis  positions.  In  the  third  and  fourth  runs, 
the  the  lower  two  dipoles  of  the  stack  were  rotated.  A  fifth 
experiment  was  also  conducted  which  considered  two  dipoles 
located  on  different  radians  (both  off  the  z-axis)  and 
oriented  in  different  directions. 

Second  Qcdei  SQUID  Gradiometsc 

The  magnetic  fields  generated  by  the  dipoles  were 
detected  with  a  SQUID  system  biult  by  the  S.H.E.  Corporation. 
The  detecting  probe  (model  BMP)  consists  of  single  weak  link 
niobium  superconducting  loop  inductively  coupled  to  a  niobium 
input  coil.  Two  niobium  screws  provide  a  superconducting 
connection  between  the  input  coil  and  the  second  order 
gradiometer  loops.  The  output  coil  is  part  of  a  R.F. 
resonant  tank  circuit.  The  BMP  probe,  as  diagrammed  in 
figure  IV-1,  is  supercooled  in  a  liquid  Helium  bath  to  4.2 
degrees  Kelvin. 
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magnetic  field  as  it  exists  on  the  upper  half  of  the 
measurement  sphere. 


Figure  III-5.  The  XY  plotted  Data.  The  data  in  the 
upper  hemisphere  is  projected  onto  the  XY-plane.  The 
resulting  contour  plot  is  an  overhead  view  of  the  de¬ 
tected  field  on  the  upper  measurement  hemisphere. 


Because  the  dipoles  were  all  located  in  the  upper 
hemisphere,  the  overhead  XY  view  provided  a  sufficient 
discription  of  the  various  field  patterns  and  the  XZ  view  was 
not  necessary.  The  IBM  control  language  used  to  run  the 
"System  II"  software  is  given  in  Appendix  B. 


The  blank  spaces  shown  in  the  file  names  above 
correspond  to  the  run#  which  the  user  enters  interactively 
when  the  program  is  executed.  The  program,  titled  MAGCAL2, 
and  its  control  language  can  be  found  in  Appendix  A. 

Plotting  the  Data 

The  data  from  both  the  simulation  runs  and  the  actual 
conducting  sphere  experiments  were  plotted  via  a  sofeware 
package  called  "Surface  II  Graphics  System"  (13).  The 
software  was  designed  to  display  spatially  distributed  data 
in  the  form  of  a  contour  map.  Although  originally  created 
for  the  purpose  of  studying  geological  formations,  the 
software  can  be  effectively  used  to  plot  any  planer 
distributed  data.  The  primary  constraints  are  that  the  data 
must  lie  on  a  plane,  be  single  valued,  and  have  locations 
described  by  two  orthogonal  coordinants. 

Because  the  Surface  II  software  requires  the  data  to  be 
located  in  a  plane,  only  two  dimensional  views  of  the 
measurement  sphere  were  plotted.  Recall  that  the  Magcal 
program  calculates  the  detected  magnetic  field  in  spherical 
coordinants  and  then  converts  to  cartesian  coordinants.  In 
order  to  project  the  spherical  data  onto  a  plane  without 
being  double  valued,  only  one  hemisphere  can  be  plotted. 
Figure  III-5  depicts  how  the  upper  hemisphere  is  projected 
onto  the  XY-plane.  Magcal  accomplishes  this  by  printing  into 
file  XYPL _  only  data  with  a  Z  component  greater  than  or 


equal  to  zero.  The  resulting  plot  is  an  overhead  view  of  the 


Read  model  parameters 


/Print  parameters  into  / 
file  SPHR _  / 


Move  to  the  next 
measurement  point  on 
the  sphere 


Convert  measurement 
point  to  cartesian  coor 


Transform  measurement 
point  to  rotated  system 
of  dipole 


Calculate  magnetic 
flux  density 


Total  field 
=  sum  of  three  dipole 
contributions 


'  Print  appropriate  coor 
and  total  flux  density 

into  files  SPHR _ ,  XYPL _ , 

and  XZPL _ 


No 

is 

— 

All  measurement 

points  evaluated? 

X  ^es 

(stoEy 

Figure  III-4.  Flowchart  of  program  MAGCAL2 
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to  the  position  of  the  first  set  of  gradiometer  coils  (recall 
equation  (10)).  It  will  be  shown  in  chapter  IV  that  this 
radius,  R,  is  equal  to  10.98  centimeters. 

The  computer  simulation  of  the  detected  field  required  a 
systematic  method  of  encompassing  the  sphere.  Naturally,  a 
spherical  coordinate  system  was  chosen.  The  equator  of  the 
sphere  was  divided  into  28  equal  units.  The  number  of 
radians  between  each  unit  is  therefore: 


2  n  radians 


28 


0.2244  radians 


(15) 


As  a  result,  a  grid  of  28  meridians  and  28  parallels  was  used 
in  the  simulation. 

The  simulation  program  consisted  of  a  Fortran  77  routine 
which  was  run  on  an  IBM  30/30  digital  computer.  The  routine 
interactively  obtained  the  simulation  parameters  peculiar  to 
each  individual  run  and  then  entered  a  loop  that  accomplished 
the  following  tasks: 

1)  Move  to  a  measurement  position  on  the  sphere. 

2)  Convert  the  measurement  point  to  cartesian  coordinants. 

3)  Calculate  the  magnetic  flux  density  by  using  the  trans¬ 
formation  of  equation  (8)  and  the  flux  density  of  equa¬ 
tion  (14)  . 

4)  Print  input  and  output  data  into  file  SPHR _ . 

5)  Print  upper  hemis.  XY  coor.  and  data  into  file  XYPL _ . 

6)  Print  lat.  "  (Y+)XZ  "  XZPL _ . 
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B 


j  j 

/  ^  /  Pa  sin  0  sin  9  / 

Jm  mm  rm  m 

m=l  n=l 


r  f  3/2(  r  ) 

n  m  n 


where , 


Magnetic  flux  density . TESLA 

Dipole  moment . AMP*METER 

Dipole  distance  from  origin . METER 

Coordinant  rotation  transformation 

R  =  origin  to  1st  set  of  coils . METER 

R+  <5  =  origin  to  2nd  set  of  coils . METER 

R+2d  =  origin  to  3rd  set  of  coils . METER 

2  =  positve  loops  in  1st  set  of  coils 

-4  =  negative  loops  in  2nd  set  of  coils 

2  =  positive  loops  in  3rd  set  of  coils 


f  (r: 

m 


1  - 


2  a  Cos  0m 

m  m 


am\  2 


Generating  iiie  Compute r  Simulation  Data 

The  pyrex  sphere  used  in  the  experiments  has  a  radius  of 
9.1  centimeters.  Note  from  equation  (14),  that  the  detected 
flux  density  is  independent  of  the  radius  of  the  conducting 
sphere.  Grynzspan's  derivations  achieve  this  result  based  on 
an  expansion  that  is  valid  only  for  the  field  outside  the 
conducting  sphere  (outside  the  source  region)  (11:22,26-29). 
Since  all  the  field  measurements,  both  simulated  and 
experimental,  were  taken  outside  the  conducting  sphere,  its 
radius  is  not  a  parameter  in  the  simulated  data. 

The  computer  generated  fields  are  represented  by  equi- 
flux  density  cotours  plotted  on  an  imaginary  measurement 


sphere.  The  radius  of  the  "measurement  sphere"  corresponds 


By  inspection,  the  first  two  components  of  equation  (12)  are 


seen  to  go  to  zero.  However,  algebraic  expansion  shows  that 
the  same  is  not  true  for  the  2nd  and  higher  order  components. 
Because  of  the  factorial  division  found  in  the  Taylor  series 
expansion,  the  higher  order  gradients  quickly  become  more  and 
more  insignificant.  Although  all  gradients  higher  than  the 
first  order  will  be  detected  by  the  gradiometer  of  figure 
III-3,  the  2nd  order  is  most  prominent  and  hence  the  name 
"second  order  gradiometer"  is  given. 

It  is  now  desireable  to  apply  the  2nd  order  gradiometer 
concepts  to  the  expanded  3-dipole  Grynszpan  model.  Since  the 
model  considers  only  a  single  scaler  component,  and  assuming 
the  model  has  similar  magnetic  characteristics  of  free  space, 
the  relation  between  the  magnetic  field  strength(H)  and  the 
magnetic  flux  density(B)  is  given  by  the  scaler  equation 
(5:260) : 


where 


B  =  u  H 

0 

u^  =  permeability  of  free  space 
=  47T  x  10_  HENRY/METER 


(13) 


Recalling  the  multipole  magnetic  field  strength  of 
equation  (9),  the  detected  flux  density  of  equation  (11),  and 
the  relation  given  in  equation  (13),  the  flux  density 
detected  by  the  2nd  order  gradiometer  is: 
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The  total  flux  detected  is  the  sum  of  the  three  flux 


density  and  the  loop  area  products: 

0  =  2B  (R)  A  -  4B  (R+<5)  A  +  2B(R+2<5)A  (10) 


where 


B(r)  =  Magnetic  flux  density  w.r.t.  position 
R  =  Radial  distance  to  first  coil 
A  =  Area  of  each  coil 

<5  =  Distance  between  each  set  of  coils 


Since  each  coil  has  the  same  area,  the  total  detected  flux 
density  is: 

B  =  2B  (R)  -  4B  (R+5)  +  2B(R+2<5)  (11) 

To  show  how  the  spacial  discriminating  gradiometer  does 
not  detect  the  constant  and  the  1st  gradient  of  equation 
(11),  perform  a  Taylor  series  expansion  about  the  point  q: 

B" (q) (R-q) 2  Bn(q)(R-q)n 

B  =  B  (q)  +  B'  (q)  (R-q)  +  -  +  ...  + 

21  n! 


B  (q) 

[  2  -  4  +  2  ] 

+  B '  (q) 

[2  (R-q)  -  4  (R+<5 -q )  +  2(R+2<5-q)J 

B"(q) 

[2  (R-q)  2  -  4{R+<5-q)2  +  2(R+2<5-q) 

2! 

B°(q) 

[  2  (R-q)  P  -  4(R+<5-q)  °  +  2(R+2<5-q)  ' 

n ! 

ICilfi  Field  Detected  a  2nd  Order  Gradiometer 

The  2nd  order  gradiometer  consists  of  a  single  wire 
wound  into  three  sets  of  coils.  The  center  set  has  twice  the 
number  of  turns  and  is  wound  in  the  opposite  sense  of  the 
other  two  sets.  Figure  III-3  illustrates  how  a  2nd  order 
gradiometer  will  detect  the  magnetic  field  at  three  locations 
corresponding  to  the  three  loops.  The  magnetic  field 
component  normal  to  the  loops  will  be  spacially  discriminated 
in  the  Oth  and  1st  order  derivatives.  In  other  words, 
constant  fields  and  linearly  changing  fields  (with  respect  to 
distance)  perpendicular  to  the  2nd  order  gradiometer  are  not 
detected.  This  can  be  shown  by  considering  the  total  flux 
detected  by  the  three  magnetometers. 


[•■TO 


Now  consider  three  dipoles  located  anywhere  in  the 


sphere.  As  long  as  each  dipole  is  tangential  to  the  sphere 
surface,  it  can  be  fit  to  a  Grynzspan  model  rotated  with 
respect  to  the  standard  cartesian  coordinant  system.  With  a 
given  measurement  point,  each  model  identifies  the  point  in 
its  own  coordinant  system  and  calculates  its  respective  field 
contribution.  The  resultant  radial  field  is  the  scaler  sum: 


H 


E 


Sin  $mSin  9m 

47r  r2  y  ’/2 


where 


(9) 
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P  = 
a  = 

7 T  = 


R  = 


T  = 


Dipole  Moment 
Tangential  dipole's 
3.1416 

Field  distance  from 
Coordinant  rotation 


distance  from  origin 
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transf ormation . 
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Note  that  the  output  of  interest  is  given  in  the  standard 
coordinant  system.  It  must  also  be  described  in  the  rotated 
system  if  the  rotated  dipole  contribution  is  to  be  calculated 
correctly.  This  is  accomplished  through  a  linear 
transformation  (21:397-398): 


Cl 

C2 

C3 

C4 

C5 

C6 

C7 

C8 

C9 

where , 

coefficients  C1-C9  are  given  by  the 
dot  product  between  the  six  axis: 
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Recall  from  equation 
in  spherical  coordinants. 


(5)  that  Grynzspan 
Transforming : 


solved  his  model 


R  =  V  x  * 2  +  ~  + 

9  =  arccos(z*/R) 

$  =  arctan(y*/  x*) 


. 
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Consider  a  second  Grynzspan  dipole  model  occupying  the 
same  sphere,  but  belonging  to  a  rotated  coordinant  system. 
This  is  illustrated  in  figure  III-2.  The  resulting  field 
generated  by  the  two  dipoles  at  some  point  P(x,y,z)  will  be 
the  vector  sum  of  H(x,y,z)  and  E{x*,y*,z*). 


z 


Figure  III-2.  A  Two  Dipole  Model.  The  Second  dipole 
belongs  to  a  rotated  Grynzspan  coordinate  system.  The 
resulting  field  will  be  H(x,y,z)  +  H  (x*,y*,z*). 


The  Pyrex  Sphere  And  Mount 

In  collecting  the  experimental  magnetic  field  data,  the 
SQUID  was  fixed.  Therefore,  the  conducting  sphere  was 
mounted  on  a  two  axis  gimbal  capable  of  systematic 
positioning. 

The  conducting  sphere  was  simply  a  3000  milliliter 
spherical  Pyrex  flask  filled  with  0.9%  saline  solution.  The 
sphere  had  a  single  openning.  The  dipoles'  mount  and  wires 
entered  through  a  #10  rubber  stopper. 

The  gimbaled  sphere  mount  was  constucted  entirely  of 
wood  and  plastic.  The  sphere  sat  firmly  in  a  rubber  lined 
wooden  bowl.  The  sphere  neck  passed  through  the  bottom  of 
the  bowl  and  through  a  tight  fitting  plastic  sleave.  The 
sphere-  bowl-sleave  assembly  moved  together  as  a  unit.  This 
assembly  was  mounted  on  a  swing.  A  short  plastic  pipe 
through  the  seat  of  the  swing  provided  a  snug  bearing  for  the 
sphere  assembly.  A  piastre  ring  with  three  nylon  set  screws 
clamped  tightly  to  the  length  of  the  neck  protruding  through 
the  bottom  of  the  swing.  Thus  the  sphere  assembly  was 
allowed  snug  rotation  but  restricted  from  backing  out  of  the 
swing . 

The  swing  assembly  was  supported  by  two  vertical 
members.  The  hinges  consisted  of  nylon  bolts  passing  through 
nylon  bearings  and  fastened  with  plastic  nuts.  The  axis  of 
swing  corresponded  to  the  horizontal  axis  of  the  Pyrex 
sphere.  The  vertical  members  supporting  the  swing  assembly 
attached  to  a  heavy  circular  base  with  a  14  inch  diameter. 


A  plastic  360  degree  protractor  was  mounted  on  each  of 
the  two  axis.  One  protractor  attached  to  the  swing  seat  and 
measured  sphere  rotation  about  the  z-axis.  A  pointer 
protruding  from  the  wooden  bowl  assembly  indicated  the 
corresponding  angle  0.  The  second  protractor  was  attached  to 
one  side  of  the  swing.  The  protractor  moved  with  the  swing 
and  measured  rotation  about  the  horizontal  axis.  A  pointer 
fixed  to  the  vertical  member  supporting  the  swing  indicated 
the  corresponding  spherical  coordinant  parameter  0.  The 
conducting  sphere  and  mount  is  shown  in  figure  IV-3. 


Figure  IV-3.  The  Conducting  Sphere  and  Mount. 
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Mount 


The  Dipole  Suppaxt 

Each  current  dipole  was  made  from  a  single  conductor 
shielded  cable  (M17/119-RG  MIL-C17E  77314  Plastoid  Corp) . 
With  the  shield  rolled  back  to  give  the  appropriate  dipole 
length,  a  small  ball  of  solder  was  applied  to  both  the  end  of 
the  single  conductor  and  the  rolled  back  shield.  In  this 
way,  current  flowed  through  the  saline  solution  to  produce  an 
effective  current  dipole. 

For  the  experiments,  the  dipoles  had  to  be  rigidly  fixed 
at  a  known  position  and  orientation.  A  plexiglass  mount 
capable  of  holding  three  dipoles  on  a  radial  stack  was  built. 
Each  dipole  could  be  rotated  independently  about  the  z-axis 
and  clamped  in  position.  Figure  IV-4  is  an  expanded  drawing 
of  the  dipole  mount.  Because  of  their  negligible  magnetic 
properties,  a  stainless  steel  screw  and  a  brass  anchor  were 
used  to  clamp  the  mount  together. 

The  stem  of  the  dipole  mount  was  placed  through  a  nylon 
disk  designed  to  fit  snugly  in  the  spherical  flask's  neck. 
Together  with  a  #10  size  rubber  stopper,  the  dipole  mount 
could  be  rigidly  positioned  within  the  conducting  sphere. 
The  mount  stem  and  the  flask  lip  were  marked  to  insure  proper 
dipole  depth  and  orientation. 

During  the  experiments,  the  SQUID  was  fixed  and  the 
conducting  sphere  was  rotated  about  the  origin.  Recall  from 
chapter  III  (see  figure  I I I — 3 )  that  the  measurement  sphere 
radius  must  be  defined  by  the  location  of  the  SQUID's  flux 
detecting  loops.  With  the  bottom  of  the  SQUID  dewar  placed 


_ Liquid 

Helium 


Figure  IV-5.  Dipole  Stack  Within  Conducting  Sphere 
and  the  SQUID  detection  Coils. 
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one  centimeter  above  the  sphere  surface,  the  bottom  detection 


coils  are  10.98  centimeters  from  the  origin.  Figure  IV-5 
shows  the  relation  between  the  SQUID  detection  coils,  the 
conducting  sphere,  and  the  dipoles. 

Experimental  Setup 

The  SQUID  sensor  hangs  from  an  aluminum  and  plexiglass 
mount  supported  by  wooden  beams.  The  entire  SQUID  support 
frame  sits  on  a  concrete  floor  in  a  corner  of  the  AFAMRL/MEG 
laboratory. 

The  support  equipment  used  in  the  experiment  include: 


Waveform  Generator  .  Wavetek  175 

Oscilloscope  .  Tektronix  7623A 

Comb  Filter  .  60  Hz,  by  SRL 

Amplifier  .  Grass  MOD.P511H 

Dual  Channel  FFT  Analyzer  .  Nicolet  660B 


A  25  hertz  sinewave  of  1.41  VAC  (RMS)  was  simutaneously  sent 
to  each  dipole  lead.  To  prevent  short  circuit  accidents  and 
to  determine  the  current  flow,  a  1000  ohm  resistor  was  put  in 
series  with  each  dipole.  The  sphere  mount  sat  on  top  a  base 
constucted  of  wood  and  cement  blocks,  and  the  SQUID  was 
positioned  one  centimeter  vertically  overhead  as  shown  in 
figure  IV-5  and  figure  IV-7. 

The  signal  detected  by  the  SQUID  passed  through  the  RF 
Head,  the  SQUID  controller,  an  amplifier,  a  60Hz  comb  filter, 
viewed  by  an  oscilloscope,  and  finally  to  channel  B  of  a 
spectrum  analyzer.  The  1.41  VAC  dipole  input  signal  provided 
a  reference  to  channel  A  of  the  spectrum  analyzer.  This 
setup  is  shown  in  figure  IV-6. 


Figure  I V-6.  The  Experimental  Setup. 
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Figure  IV-8.  The  30X  SQUID  Controller. 


Figure  IV-11.  Tektronix  7623A  Oscilloscope. 


Figure  IV-12.  Nicolet  660B  Dual  Channel  FFT  Analyze 


the  Data 


In  the  first  four  experiments,  the  dipole  stack  mount 
was  set  and  clamped  to  the  appropiate  orientations.  The 
fifth  experiment  did  not  consist  of  a  dipole  stack  and  hence 
a  separate  means  of  supporting  the  dipoles  was  used.  Figures 
IV-13  and  IV-14  show  the  dipole  mount  for  the  first  four  and 
the  fifth  experiments  respectively.  Figures  IV-15  through 
IV-19  illustrate  the  five  multipole  models  examined. 


Figure  IV-13.  The  Dipole  Mount  used  in  Experiments  #1-4. 
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Figure  IV-14.  The  Dipole  Mount  used  in  Experiment  #5. 


Figure  IV-15.  Experiment  #1  Model. 
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Figure  IV-16.  Experiment  #2  Model. 
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Figure  IV-17.  Experiment  #  3  Model. 
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Figure  IV-18.  Experiment  #4  Model 
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Figure  IV-19.  Experiment  #5  Model 


Figure  V-9.  Calculated  Field  for  Model  #5. 
The  measurement  sphere  radius  is  10.98  cm. 
Magnitudes  are  given  in  Tesla  XEE-11. 
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Figure  V-8.  Experiment  Field  for  Model  #4. 
The  measurement  sphere  radius  is  10.98  cm. 
Magnitudes  are  given  in  Tesla  XEE-11. 
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Figure  V-7.  Calculated  Field  for  Model  #4. 
The  measurement  sphere  radius  is  10.98  cm. 
Magnitudes  are  given  in  Tesla  XEE-11. 
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Figure  V-6.  Experiment  Field  for  Model  #3. 
The  measurement  sphere  radius  is  10.98  cm. 
Magnitudes  are  given  in  Tesla  XEE-11. 
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Figure  V-4.  Experiment  Field  for  Model  #2. 
The  measurement  sphere  radius  is  10.98  cm. 
Magnitudes  are  given  in  Tesla  XEE-11. 


67 


66 


Figure  V-2.  Experiment  Field  for  Model  #1. 
xue  measurement  sphere  radius  is  10.98  cm. 
Magnitudes  are  given  in  Tesla  XEE-11. 


Figure  V-l.  Calculated  Field  for  Model  #1 
The  measurement  sphere  radius  is  10.98  cm. 
Magnitudes  are  given  in  Tesla  XEE-11. 


When  model  #5  was  evaluated  in  the  conducting  sphere 
experiment,  the  "primary  null  +  branches"  field  did  not 
completely  form.  Note  from  figure  V-10  that  the  positive 
flux  peak  belonging  to  the  -y  oriented  dipole  did  not  appear. 
However,  the  contour  lines  in  the  vacinity  do  indicate  the 
peak's  presence.  By  inspecting  the  calculated  field  of 
figure  V-9,  it  is  apparent  that  this  missing  peak  should 
have  appeared  at  the  coordinant  (-5,4).  It  is  believed  that 
imperfections  in  the  experimental  model  plus  sphere 
inhomogenetxes  resulted  in  the  positive  peak  being  swallowed 
up  by  the  neighboring  negative  peaks. 

The  data  used  to  plot  the  computer  generated  fields  and 
the  actual  fields  measured  with  the  second  order  SQUID 
gradiometer  are  given  in  appendixes  C  and  D  respectively. 
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the  fields  took  on  a  twisted  appearance  as  compared  to  the 
first  two  models.  The  null  line  favors  the  shallow  dipole 
at  field  points  directly  above  the  dipole  stack,  and  favors 
the  deeper  dipoles  at  field  points  away  from  the  dipole 
stack. 

The  fields  plotted  for  the  first  four  experiment  models 
are  highly  distorted.  The  primary  cause  for  this  distortion 
is  believed  to  be  related  to  the  dipole  mount.  The  3rd  and 
4th  experiments  show  more  distortion  since  the  rotated  mounts 
interfere  to  a  higher  degree  with  the  current  densities.  The 
dipoles  in  the  5th  experiment  were  held  in  place  by  a  far 
less  bulky  mount  and  thus  caused  the  least  distortion. 

The  5th  model  considers  the  field  generated  by  two 
dipoles  located  on  different  radials.  The  calculated  field 
was  found  to  have  a  branched  null  line.  The  primary  null 
followed  the  average  direction  of  the  contributing  dipoles 
with  the  branches  extending  from  the  primary  null  and 
passing  over  the  vacinity  of  respective  dipoles.  See  figure 
V-9.  An  additional  computer  generated  plot  (not  included  in 
this  report)  showed  that  as  the  dipoles  of  model  #5  were 
moved  closer  to  the  sphere  origin,  the  null  branches  moved 
closer  to  the  primary  null.  For  deep  model  #5  dipoles,  the 
branches  fused  with  the  primary  null,  thus  eliminating  the 
two  close  proximity  density  peaks  and  resulting  in  a  plot 
resembling  that  of  a  dipole  stack.  This  result  is  expected 
since  the  dipoles  will  meet  in  the  limit  of  reaching  the 
origin. 
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V.  Results 


Recall  the  five  models  shown  in  figures  VI-15  through 
VI-19.  Two  plots  were  drawn  for  each  of  the  five  models: 

1.  Calculated  field 

2.  Experiment  field 

Each  plot  is  an  overhead  view  of  the  equi-magnitude 
magnetic  flux  density  contours  that  encompass  the  measurement 
sphere.  An  overhead  diagram  with  the  parameters  describing 
the  various  dipole  locations  and  orientations  is  included 
with  the  calculated  field  plot.  See  figures  V-l  through  V-10. 

When  comparing  the  computer  generated  field  with  the 
conducting  sphere  experiment  field,  both  similarities  and 
discrepancies  are  noticed.  In  each  model,  first  note  the 
similarity  of  the  general  form  of  the  two  fields.  A  null 
contour  separates  the  positve  flux  densities  from  the 
negative  flux  densities  and  maximum  density  peaks  appear 
within  an  expected  region. 

The  data  of  the  experiment  plots  are  on  the  same  order 
of  magnitude  but  about  three  times  larger  than  the  calculated 
data.  This  result  indicates  that  the  measurement  sphere 
radius  used  in  the  calculated  fields  may  have  been  too  small. 
This  is  a  real  possibility  since  the  actual  radius  was 
difficult  to  measure  and  the  magnetic  field  varies  as  the 
square  of  the  distance.  Since  this  factor  remains  consistent 
throughout  the  data,  important  contour  shape  information  is 
preserved. 

As  the  dipoles  were  rotated  in  the  3rd  and  4th  models, 


The  polarity  of  the  detected  flux  with  respect  to  the 
input  voltage  was  determined  by  the  phase  shift  between  the 
dipole  input  voltage  and  the  SQUID  output  voltage.  Recall 
from  chapter  II  that  the  output  voltage  lags  the  input 
voltage  for  positive  flux  (w.r.t.  input)  and  leads  for 
negative  flux  (w.r.t.  input).  For  each  measurement  point, 
the  660B  also  calculated  the  transfer  fuction  of  the 
output/input.  In  the  collection  of  the  data,  corresponding 
phase  shifts  where  primarily  clustered  into  two  groups 
approximately  180  degrees  apart.  Phase  shifts  that  deviated 
from  the  rule  greatly  were  associated  with  very  low  signals 
and  were  assumed  to  be  near  the  null  contour.  Using  the  right 
hand  rule  as  a  guide,  positive  phases  were  assigned  negative 
flux  densities  and  negative  phases  were  as signed  positive 


In  each  of  the  five  experiments,  the  conducting  sphere 
was  placed  one  centimeter  directly  beneath  the  SQUID 
dewar.  With  both  sphere  orientation  angles  set  to  zero,  the 
linear  axis  of  the  SQUID  was  set  to  coincide  with  the  z-axis 
of  the  sphere. 

Data  collection  was  accomplished  by  systematically 
rotating  the  gimbal  mounted  sphere  to  spherical  grid  points 
defined  by  the  horizontal  and  vertical  protractors.  The  grid 
consisted  of  163  points  in  the  northern  hemisphere  only.  The 
measurement  points  were  converted  to  cartesian  coordinants 
and  their  resulting  x-y  values  are  given  in  the  experimental 
data  listings  of  appendix  D. 

For  each  point,  the  660B  spectrum  analyzer  averaged  5 
output  voltage  measurements  over  a  span  of  20  seconds.  The 
following  factors  were  used  to  convert  the  detected  RMS 
voltage  to  magnetic  flux  density  (Ref  22) : 

SQUID  sensitivity  on  xlO  scale  ..  1.74  xlO  7gaus/volt 

Grass  amplifier  .  xlOO,  +134  shift  @  25  Hz 

60  Hz  Comb  filter  . .  xl ,  +2  shift  @  25  Hz 


where , 

v  t  =  output  of  330X  control  unit 
v2  =  output  of  Grass  amplifier 


Figure  V-10.  Experiment  Field  for  Model  #5. 
The  measurement  sphere  radius  is  10.98  cm. 
Magnitudes  are  given  in  Tesla  XEE-11. 
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VI.  Discussion  and  Conclusions 


The  model  developed  in  this  study  considers  the  simple 
case  of  a  homogeneous  conducting  sphere.  The  external  radial 
magnetic  field  component  existing  outside  the  sphere  results 
from  the  dipole  currents  alone.  Cuffin  and  Cohen  (12:374- 
375,  381)  showed  that  the  volume  currents  have  negligible  or 
zero  contribution. 

The  perfect  homogeneous  sphere,  however,  could  not  be 
replicated  in  the  actual  conducting  sphere  experiments.  The 
presence  of  the  dipole  wires,  mount,  etc,  does  not  allow  the 
volume  currents  to  symmetrically  cancel.  This  was  most 
prominent  in  the  cases  where  a  large  mass  of  nonconducting 
material  was  present  near  the  active  dipoles.  As  a  result, 
the  unsymmetric  volume  currents  did  contribute  substantially 
to  the  detected  field  and  warped  it  with  respect  to  the 
homogeneous  model.  Here,  it  should  be  noted  that  the  degree 
of  field  warping  is  dependent  on  the  proximity  of  gross 
inhomogeneties  to  the  active  dipoles. 

HuLtipole  Ariay  Zisld  Estimator 

Recall  that  the  desired  goal  of  modeling  the  magnetic 
field  activity  of  the  cerebral  cortex  is  to  find  the  simplest 
description  of  the  sources  that  produce  the  observed  field 
pattern.  This  may  be  achieved  by  modeling  the  cortex  as  an 
array  of  stacks  containing  one,  two,  or  three  dipoles  each. 
Cortical  folds  rich  in  pyramidal  cells  tangential  to  the  head 
surface  would  be  expected  to  have  a  larger  stack  than  surface 


cortex  and  thus  have  a  greater  contribution  to  the  detected 
field. 

The  parameters  required  to  describe  each  dipole  must 
include  dipole  moment,  distance  from  the  origin,  and  a  3X3 
transformation  matrix  defining  dipole  orientation.  With  11 
parameters  required  for  each  dipole,  the  number  of  model 
variables  can  grow  to  an  enormous  size.  Since  this  size  can 
be  reduced  when  modeling  isolated  sources  by  setting  most  of 
the  dipoles  to  zero,  the  overall  cerebral  cortex  model  may  be 
developed  from  the  summed  effects  of  the  reduced  models. 

Of  course  the  complete  functional  computations  of  the 
cortex  are  far  too  complex  to  consider  presently.  It  is, 
however,  reasonable  to  develop  models  of  observed  fields 
with  the  dipole  parameters  describing  isolated  brain  states. 
If  a  given  state  can  be  modeled  reliably,  it  may  be  possible 
to  apply  stochastic  estimation  techniques  successfully. 

For  example,  consider  the  estimation  of  2  states  within 
the  motor  cortex  corresponding  to  the  flexion  and  extension 
of  the  index  finger.  Each  state  must  be  described  by  several 
parameters  indicating  all  pertinent  dipole  contributions. 
Although  isolated  state  fields  may  be  obtained  through 
repetitive  averaging  experiments,  specific  dipole  locations 
are  difficult  to  determine  in  the  case  of  a  multipolar 
source.  State  dynamics  and  noise  characteristics  also  need 
to  be  fully  understood. 

The  uncertain  model  parameters  may  be  estimated  through 
maximum  likelihood  techniques.  The  likelihood  is  a  function 
of  the  parameters  to  be  estimated  and  the  measured  data.  The 


goal  is  to  find  the  value  of  the  parameter  vector  that 


maximizes  the  likelihood  function.  These  techniques  are 
covered  in  detail  in  chapter  10  of  reference  23. 

At  this  point  it  is  not  even  known  if  state  and 
parameter  estimation  is  possible.  Many  experiments  need  to 
be  accomplished  before  the  models  may  even  be  considered  an 
adaquate  representation. 

Recommendations  for  Continuing  Studies 

Data  collection  using  a  single  second  order  SQUID 
requires  hours  of  work  before  a  field  plot  can  be  generated. 
If  just  one  SQUID  is  available,  it  is  suggested  to  automate 
data  collection  with  a  digital  computer.  This  could  be  done 
by  using  the  output  ports  of  the  660B  spectrum  analyzer  to 
transmit  magnitude  and  phase  data  directly  to  a  plot  file. 
Although  the  conducting  sphere,  or  subject,  would  have  to  be 
manually  placed,  the  plot  turnaround  time  could  be 
significantly  reduced. 

The  most  desirable  system  would  allow  all  data  points  to 
be  taken  simutaneously .  A  multi-array  SQUID  would  surely 
enhance  these  experiments.  In  the  interest  of  studying  real 
time  cerebral  fields,  the  desired  magneto-encephalography 
analyzer  should  include  an  array  of  at  least  100 
semiconductor  SQUIDs. 

Conclusion 

This  project  took  Grynszpan's  dipole  in  a  sphere  model 
and  expanded  it  to  include  any  combination  of  three 


tangential  dipoles.  The  results  show  that  dipoles  oriented 
on  a  stack  do  not  produce  multiple  peaks,  but  rather  a 
twisted  version  of  a  single  dipole  field.  Dipoles  located  on 
different  radials  will  produce  multiple  peaks  provided  the 
dipoles  are  separated  by  a  measurable  distance. 

The  experiments  using  an  actual  conducting  sphere 
support  the  model  results.  The  pyrex  mount  used  to  support 
the  dipoles  caused  a  disturbance  in  the  volume  currents.  It 
was  found  that  the  greatest  field  distortion  occured  when 
large  masses  of  the  nonconducting  pyrex  were  near  the  active 
dipoles . 

Mathematical  modeling  of  cerebral  cortex  activity  is 
still  in  its  infancy.  The  rapid  advancement  of  multiarray 
SQUID  technology  is  necessary  to  develop  experiments  needed 
^  *  to  determine  model  adaquacy. 
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Appendix  A 


This  section  contains  the  Fortran  77  routine  and  the  IBM 


30/30  control  language  used  to  generate  the  calculated  field 
data.  Also  included  is  an  example  execution  of  the  program. 


Job  Control  Language  for  MAG  CAL  2 


000 10  PROC  1  SERIES 

00020  CF REE  ATTR(XXX)  FI(FT1 1F00I ,  FT  1 2F00 1 , FT ! 4F00 1  ) 

00030  ATTR  XXX  RECFM(F  B)  LRECL ( 80 )  BLKSIZE( 12960 ) 

OCO^O  ALLOC  FI(FT1 1F001 )  DAtSPHRASERIES . )  NEW  TRACK  SF(2,2)  USING(XXX) 

000S0  ALLOC  FKFT12F001)  DA ( XYPL&SERIES . )  MEW  TRACK  SP(2,2)  USING(XXX) 

00060  ALLOC  FI(FTHFOCI)  DA(XZPL4SERIES. )  NEW  TRACK  SP(2,2)  USING(XXX/ 

00070  CALL  (MAGCAL2) 

END  CF  DATA 


-CLIST  Data  ifl  Control  MAG  CAL  2 


00010  //HETPJD2  JOB  HET , DEREGO , NOTIFY'HETPJD , TIME' ( 2 ) , MSGCLASS'Z 

00020  //  EXEC  FORTOCL 

00030  //SYSPRINT  DD  SYSOUT'Z 

00040  //FORT .SYSIN  DD  DSN=HETPJD . FORT ( MAGCAL2 ) , DISP'SHR 

000S0  / /LKED .SYSLMOD  DD  DSN'HETP JD .L0AD(MAGCAL2 ) , DISP= ( NEW , CATLG ) , 

00060  //  UNIT-DISK, VOL-SER-PACK01,SPACE-(TRK, (10, 1G, 2)), 

00070  //  DCB'SYSI .LINKLIB 
END  OF  DATA 


Broqram  magcal2 


our.;:':  . 

.  i  o  ij ;  ii  c 
onrur.  r 
,nu.;'.n  c 

OJiUO  C 
,.iilU7n  C 
n nr, so  c 
000311  c: 
oui  nr.  c 
0011 0  c 

00120  C 

00  I  00  c 
onion  r: 
00150  C 
00160  C 
00170  C 
001 130  C 
00130  C 
002 or.  c 
00210  C 
00220  C 
00230 
00240  C 
00050  C 
00250  C 
00270  C 
JO 230  C 
00290  C 
00 3 00  C 

00310  C 
On  3 2 ! l  !.' 
00330  3 
on  :;40  r. 

JO  060 
00350  C 
00  3  70  .. 
Ol:  3211  r 
no .  'in  i. 

,'i  i.i  <41 1  i.i  r. 

IIUll'l  I' 
i -0  1 1  t; 

.T'-i  ;o  r 

0;, 4 4i.i  c 
..!!’>  *lc.li  3 


ilUL  •’Iylf'Oi.  r.  I  ."I  i  ;.PrlErCE 

iuT  PAUL  0.  DE  REDO 
5F  I  T  .  UPl.fB.  OHIO 
OCTOBER  30.  1 jS4 

THE  PROGRAM  CALCULATES  THE  MAGNE1IC  KlUX 
DENSITY.  AT  A  UNIFORM  RADIUS.  DETECTED  BY  THE  AFMRL 
SQUID  GRAD  IOMETER  OUTSIDE  A  CONDUCTING  SPHERE  CON¬ 
TAINING  3  CURRENT  DIPOLES. 

THE  CALCULATIONS  USED  ARE  BASED  ON  GRYNSZPAN ’ S 
RADIAL  FIELD  EQUATION.  REVISIONS  INCLUDE  AN  EXPAN¬ 
SION  TO  INCLUDE  THREE  TANGENTIAL  DIPOLES  POSITIONED 
ANYWHERE  IN  THE  CONDUCTING  SPHERE. 


CONCEPT IONALLY .  THE  GRYNSZPAN  MODEL  IS  ROTATED. 
COMPUTATIONALLY.  THIS  IS  DOME  BY  DESCRIBING  THE 
MEASUREMENT  POINT  OF  THE  STANDARD  CARTESIAN  COORDINATE 
SYSTEM  IN  THE  ROTATED  CARTESIAN  COORDINATE  SYSTEM. 

THE  POINT  IS  TRANSFORMED  TO  SPHERICAL  COORDINATES 
W.R.T.  1  HE  ROTATED  SYSTEM  AND  THEM  THE  FIELD  IS  CAL¬ 
CULATED  USING  GRYNSZPAN ’ S  EQUATION. 

TIE  MEASURING  POINT  IS  TRANSFORMED  FOR  EACH  DIPOLE 
SEPARATELY.  THE  RESULTS  ARE  SUMMED  AND  STORED  AS  THE 
EIElD  OF  THE  POINT  W.R.T.  THE  STANDARD  SYSTEM  IN  BOTH 
SPHERICAL  AND  CARTESIAN  COORDINATES. 

THE  INPUT  AMD  OUTPUT  Data  ARE  RECORDED  IN  A 
SINGLE  FILE.  FROM  THE  DATA  TWO  ADDITIONAL  FILES. 

FIELD  ABOVE  XY  PLANE  AND  FIELD  ABOVE  XI  PlANE  . 

ARE  CREATED.  THESE  TWO  FILES  ARE  COMPATIBLE 
WITH  THE  SURFACE  II  CONTOUR  PROGRAM. 


VARIABLE  KEY 
INPUT: 

P  =  DIPOLE  MOMENT  (AMP -Ml 

A  -  DISTANCE  BETWEEN  TANGENTIAL  DIPO-E  ORIGIN  <M) 

R  -  MEASUREMENT  SPHERE  RADIUS  (1ST  GRAD  COILS)  (M) 
C1-C3  .  3X3  TRANSFORMATION  MATRIX 

OUTPUT: 

1.  MEASUREMENT  SPHERE  FILE 

2.  XY  PLANE  PROJECTION  FILE 
2.  XZ  FVAUE  PROJECTION  FILE 
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X  ,  Y  .  i  .  .  y.  AR  :  COOK  PARmMl  Tr.Ro . 1  *.M) 


INTEGER  T.F.I 

REAL  P f  3 ) . A ;  3 ' . R . C 1 1 3 ) . C2 ( 3 ) . C3 (3) . Cd <3» . C5 < 3 • . C6 (  3) . 
407(3) ,CS ( 3 > .03(5 ) 


DO  20  1-1.3 

PRINT  •‘.’ENTER  DIPOLE  *'. I ,' MOMENT  :  X.XXE  XX  (AMP-*M)' 
READ  (B.IO)P(I) 

FORMAT (E8 .2 ) 

CONTINUE 

DO  30  1*1.3 


00670 
00680 
00690  60 
00700  30 
00710  C 
00720 
00730 
00  740 
00  750 
00  760 
00770 
00730 
00  790 
0  0300 
0  0810 
00820 
00830 
00340 
00350 
00860 
.  1 0  3  7  0 
00380  70 


PRINT  ‘.'ENTER  DIPOLE  . I . VQCN:  X.X  (M  FROM  ORIGIN)' 
READ ( 5.60 )A( I ) 

FORMAT ! F6 . 4  ) 

CONTINUE 

DO  "0  1=1.3 

PRINT  ‘.'TRANSFORMATION  MATRIX  FOR  DIPOLE* ’.I 
PRINT  ■ 

PRINT  *. '  C!  C2  C 3  ' 

PRINT  C4  C5  C6  ' 

PRINT  •».'  C7  C3  C9  ' 

PRINT  *.'  ' 

PRINT  ‘.'ENTER  THE  FIRST  ROW.  Cl  C2  C3.  IN  THIS  FORMAT:’ 
PRINT  *- .  ’  X.XXX  X.XXX  X.XXX' 

PE  AD  ( 5 .  /  0  )  1  <  1 1  ,  C2 1 1 )  .  CO  ( I ) 

PRINT  *. 'ENTER  THE  SECOND  ROW.  C4  C5  Cb .  IN  THIS  FORMAT:' 
PRINT  >.'  X.XXX  X.XXX  X.XXX' 

READ ( 5 . 70 )C4( I ) . CS ( I !  . C6 ( I ) 

PRINT  ..'ENTER  THE  THIRD  ROW.  C7  CS  CS.  IN  THIS  FORMAT." 
PRINT  «. .  •  X.XXX  X.XXX  X.XXX' 

READt  5 . 70 ' C7 ( D , C3 ( D  .CSiI) 

FORMAT (  It  7.5) 


.11.13':' I! 

o 

:i  o  a  i  n 
0 1!'j3|! 

U  3 :  •  i 

,;i  o  9  a  o 

.'i  0  s  £.  i'i 

0  n%i.. 

00370 
00330 
00390 
01000 
01010 
0102C 
0'030 
01040 
01050 
01060 
01070 
01080 
01090 
01  100 

0  til  0 
01120 
on  30 
01140 
0  1-1 50 
01160 
01 -170 
01160 
01  190 
01200 
01210 
01220 
.1 1230 
01240 
01250 
01260 
012  7  0 
01280 
0  1230 
01  H'1- 


40  CON  T  I. ‘II  IE 

PRINT  ‘.’ENTER  THE  NCOS  RADIUS  l  1ST  OR  AD  .'Oil. ) :  2 . '!  CD' 
RFf-.D  i  6  1  E 


CPE  fit  12' 

OPEN'  I  4 ) 

WRITE' 1 1 . 1 1 12) 

WRITE) 1 1 .2010) 

WRITE',  i  1  ,2020  )P(  1  )  ,P(2)  .P(3) 

WRITEt 1 1 .2030 !A( 1 ) . A < 2  > ,6(3) 

WRITE' 1 1 . 11 12) 

WRITEt I  1 .2040 ) Cl ( 1 ) . C2 (!) , C3 (I ) . C 1 ( 2 ) . C2 ( 2 ) . C3 ( 2 ) . 
-'•Cl  ( 3 ) . C2 ( 3 ) . C3 ( 3 ) 

WRITE ( 1  1 . 204 1 )C4 ( 1 ) . C5 ( 1 ) .C6( 1 ) .C4( 2 ) . C5t 2) . C6(2 ) , 
+C4f3) .CS(3) .06(3) 

WRITEt 11 .2042 ) C7  t 1 ) . C8 ( 1 ) . C9 <  I  ) . C7 ( 2 ) . C8 ( 2 ) , C3 ( 2 ) . 
+C7(3> . C3  f  3 ) . C3( 3 ) 

WRITEt 1 1 .1112) 


mo 

2010 

2020 

2030 

2040 

2041 

2042 
2050 
1  1  1  1 


1  1 


WRI Tc (  1  1 
WRITEt 1 1 
WRITEt  I  1 
WRITEt 1 1 
WRITEt  I  1 
WRITEt 1 1 
FORMA Tf ’ 


2050 )  R 
1112) 
1110) 
1112) 
111!) 
Ill, 


•HER  DATA  POINTS:  THETA=T  *0 . 2244  :  PHI«r-*Q. 


244 


FORMAT  (11a.  ' DIPuLE#1  Ml 


'DIPOLE#, 


'DIPOLE#':;' ) 


FORMAT! ’MOMENT 
FORMAT! 'LOCATN 
FORMAT' ' TRAMS 
FORMAT ( 'MATR IX 
0  R  MAT  ( 

FORMAT ('THE  ME AS 
FORMAT i 4X ,  '  BD (  T 
'  'FLUX  DENS. 

FORMA T( ' ’ ) 


,E8. 
,F7 . 


2. 12X, 


1 ':) 


EX 


E8.2. 12X.E3.2.3X 
F7 .3.1.3X.F7.3.6X 
3FS. 3.2X.3F6.3) 

3F6 . 3 . 2X . 3F6 . 3 i 
3FS . 3 . 2X . 3F  6 . T  > 

P HERE  RAD IUS  TS ' . 2X . Fb . 4 . 2X . ' M 
R )  '  .  1  OX .  '  X ' . 7X . ' Y ' . 7X . ' 7 ' ,SX. 
XEE-11  TESLA') 


(  AMP  *M ) ' ) 


(M  FR  ORN)’) 


.  3F6 .  .3 . 
.  3FF. .  3  , 


^  A  . 

'"I  \/ 


C 


“UPPER  POLE*  “ 


■  ii  ?::o 
j  1  3  4  0 
.J  i  5 1  i  I. 

•  i i  :r, n  ;■ 
:n370 
0 1  350 
0  1  330 
01400  11. 
01410  II. 
0 1 420  C 
01430  C 
01440 
01450 
01460 
01470 
01  48C 
01490 
01500  C  • 
01510 
01520 
01530  C 
01540  C 

01550 
01560 
01570 
015S0 
01590 
01-600  C 
01610  C 
01620 
01630  C 
01640  C 
01650 
01660  C 
01670 
01630 
0  I  690 
01700  C 
017!  0  C 
01720 

0  1 730 
01  740 
017511  r: 


.  A .  L  f  I E D  !  C  .  F  .  R  ,  9  .  R  .  C  I  .1  2 .  0  3  .  64  .  C:. .  6E> .  C  7  ,  CO  .  C3 ) 

» ■*■  > h'r'K1  OL »_  L u AND  11i_ R  I L)  j.Ai  15  *  *  * 

[:  j  !  :  j  1 

DO  ’  12r  f'-n.j.- 

CliL  .  ■  !ELD<T..*  . R  .  A.  P  .11  .  02 . 63 .  C4 . 06  .  CS .  07  . C3.C9  ) 
'0  CONVINCE 
1  CONTINUE 

•*  *■  ■  * i_ 0 14 E . \  POLE4-*-* 

F-D 
T  =  1  4 

CALL  FIELD! T .F.R.A.P. Cl .C2.C3.C4.C5.C6.C7.CS.C9) 
CLOSE (  1 1  ) 

CLOSE  M2) 

CLOSE < 14) 


SUBROUTINE  FIELD< T.F .R.A.P.C1 .C2.C3.C4.CS.C6 . C7 . C3 . C9 ) 
INTEGER  T.F.M.N.I 

REAL  R .  A  (  3 )  .  ?  <  3 ) .  B .  D  .  RD  (  3 )  .  K  (-3  )  .  FR  <  3 )  .  SUM  1  .  SUM2  . 

+  SI (3) ,S2(3) .BD.X.T . Z.C1 (3) .C2(3> . C3( 3  > .C413) ,C5(3) .C6(3) 
+  C7 ( 3 ) . C3 ( 3 ) . C9 ( 3 ) . XX . YY . Z7 . T  T . FF 

RADIANS  BETWEEN  ALl  PARALLELS  AND  MERIDIANS 
B* .2244 

DISTANCE  BETWEEN  GRADIOMETER  COILS  (METERS) 

D  =■  .03175 


:  i  -  <  -  4 ) 


TRANSFORM  MEASUREMENT  POINT  TO  CARTESIAN 
X  -  R  *  COS ! r*  B  '  *  SII.1  T*-E. 1 
'<  =  n  >  s  i ;  j  '  p  *  Ei  >  -*  sni'T'B) 

Z  =  R  •-  CO'.-,  i  T*[:i 

DO  1127  M- 1  . 


I  !,'■/!  I 


,1 !  'I'll 

’ KAfL r  ora  •“  .v  o:  r'Oini  !>;j  mr  o  osiem 

,j ; ui 

r:  :  i  ;i  •  *  ;  *  .  j  i  ii  >  *  y  *  :  ( ;i  >/ 

i.i  i  :vim. 

> v  --  r-,4( i  f  i m  >  <  v  i  r.G(h'<2 

.1 ’ Ml 

22  -  .  MM)  >:  +  ,:m:i)*y  +  :"9t!H>Z 

(J  i 

,1 ! M  l 

. 

’  HA  ill  r  i.'KM  i F.  IF  H  I Ai_  >_OOR  U.R.T.  ROTATED  jYoIc.il 

0 1  MM 

11  -  ALUM 

u  i  8~n 

!i  nc(.:u.u  .AND.  YY.Ea.il)  THEN 

018  SO 

FT  =  n 

j 1  o  7 0 

EjSE 

01380 

ff  -  ATAN2M  yy  i 

j  !  ,3  9  ij 

END  IF 

mono 

Y 

.i  1 3i  a 

< ' 

THE  '.‘.2 C3ND  ORDER  ORAL)  I  OilL  TER  EFFECT 

oio;?u 

DO  1 1  2 6  r i  =  i  .  O' 

0 1  330 

RDi  M  i  -  R  ■■  L» » 1  .'i-  f  ) 

01000 

FR(N>  *  1  -  2*-AiN)m:05(  TTi/RD(Ni  +  ( A  (  N  )  /RD  i  N )  )  *-‘-2 

01 900 

9  1  (  N )  =  FIN)/  v  R  D  (  N  !  2  -*FR(N  )-»->!  .5) 

01060 

1 1  28 

CONTINUE 

01970 

SUM1  =  SIM)  +  5 1(2)  +  SI  (3) 

01930 

C 

0  1990 

5  2  ( M  i  =  PCI)  -<■  AIM)  -<•  SINiTT)  •*  SIN(FF)  ->  SUM1 

02000 

1127 

CONI  III UE 

02010 

BD  »  S2 M )  +  S2'2l  +  S2(3> 

02020 

c 

02030 

RECALL  THAT  THE  CALCULATED  RESULT.  BD.  HAS  EE-7  TACKED 

02030 

c 

ONTO  IT.  THIS  FACTOR  WAS  CARRIED  THROUGH  FROM  THE 

02050 

c 

PERMEABILITY  CONSTANT. 

02060 

c 

02070 

BD  =  BD  10000 

02080 

c 

TO  MAKE  PLOTTING  LESS  DIFFICULT.  BD  IS  SCALED  ANOTHER  4 

02090 

ORDERS  OF  MAGNITUDE.  THE  SCALING  FACTOR  IS  ADJUSTED  TO 

021  CO 

c 

EE-11.  SINCE  THIS  FACTOR  IS  ALSO  THE  SAME  ORDER  OF  MAG 

02110 

8 

AS  THE  EXPERIMENT  SENSITIVITY.  DATA  COMPARISON  IS  EASIER 

02120 

c. 

02130 

L 

THE  CARTESIAN  COORDINATES  ARE  PLOTTED  IN  CENTIMETERS: 

02140 

X  «  x  *  ICO 

02150 

Y  =  Y  ♦  100 

02160 

Z  =  Z  *-  100 

02170 

c 

02  1  30 

c 

spherical  data 

,j2  r  90 

WRITE!  M.l  129) .T.F.K.X.Y.Z.BD 

02200 

1 1 29 

FORMA  1  ( 2X .  ’  BD  (  '  ,  IE  .  ’  .  ’  .  12 .  ’  . '  .  Ft. .  4 .  ’  1  ' .  3F3 . 3 . 4X .  E 1 S .  S  i 

.122  1  0 
.12220 

( 

XY  PLANE  PROJECTION 

0  :■  2  :  i'i 

IF(Z.GE.O)  THEN 

022  4  C 

WRITE  ME.  11  30  .X.Y.BD 

02250 

END  IF 

02260 

C 

02270 

c 

XZ  PLANE  PROJECTION 

02230 

I’D 

02290 

IF(Y.GE.O)  THEN 

02300 

1  =  1 

02310 

END  IF 

02320 

IF(Z.GE.O)  THEN 

02330 

I  =  IU 

02340 

END  IF 

02350 

IF(  I . E 0 . 2 >  THEN 

02360 

WRITE!  14. 1 130  .X.Z.BD 

02370 

END  IF 

02380 

L 

02390 

1  130 

FORMAT ( 2F8 . 3 . 5X . E 1 5 . 3 ) 

02400 

c 

02410 

END 

END.  OF 

DATA 

83 


Example  Execution  .of.  MAG CAL 2 

The  job  control  language  JCL (MAGCAL2 ) ,  control  file 
CLIST (MAGCAL2 ) ,  and  the  routine  FORT (MAGCAL2 )  are  compiled 
by  submitting  JCL (MAGCAL2 ) .  Once  the  program  has  been 
compiled,  it  can  be  executed.  The  program  will  prompt  the 
user  for  input  data.  The  position  parameter  series  indicates 
the  particular  run  #.  For  example,  the  parameter  series  05 
will  create  and  store  the  calculated  data  in  files  SPHR05, 
XYPL05,  and  XZPL05 . 


PL. 


ENTER  POSH  U 
.16 

ENTER  DIPOLE 

iNtti. 

pa.ramftef 

:.ER  Ico 

4 1 

1 

MOMENT  : 

v  vvc  yv 

A  •  /\AL  a  A 

(  AMP-H1 ) 

I . 10  E-0S 

ENTER  DIPOLE 

4 

"> 

MOMENT: 

X.XXE  XX 

(AMP*M) 

1.10  E-OS 

ENTER  DIPOLE 
U 

ENTER  DIPOLE 

if 

3 

MOMENT: 

\/  \/\/r~  \/x/ 
A  «  A  a  CL  A  /\ 

( AMP»M) 

-  I 

4 1 

1 

L0CN:  X 

X  ( M  FROM 

ORIGIN) 

0 .070 

ENTER  DIPOLE 

if 

2 

L0CM:  X 

X  (M  FROM 

ORIGIN) 

0 . 0  7  0 

ENTER  DIPOLE 

if 

o 

L0CM :  X 

X  Cl  FROM 

ORIGIN) 

0 

TRANSFORMATION  MATRIX  FOR  DIPOLE#  I 

Cl  C2  C3 
C4  CS  CS 
C7  C3  C 9 

ENTER  THE  FIRST  ROW.  Cl  C2  C3.  IN  THIS  FORMAT: 
X.XXX  X.XXX  X.XXX 
H.906  0.000  -0.42:: 

ENTER  THE  SECOND  ROW,  C4  CS  CS.  IN  THIS  FORMAT: 
X.XXX  X.XXX  X.XXX 
0.C00  1  .00  0  C.OOt. 

ENTER  THE  THIRD  ROW.  C7  C3  C9,  IN  THIS  FORMAT: 
X.XXX  X.XXX  X.XXX 
0.420  0.000  0.906 
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TRANSFORMATION  MATRIX  FOR  DIPQLt:* 


•Cl 

C2  C  3 

CL 

CS  C6 

C7 

Co  C9 

ENTER  1 

HE  FIRST 

ROW. 

X.XXX 

X  .  X  X )  * 

\.  \/\/\ 
A  •  /•  /\  /'. 

0 . 000 

-0.9  of,  - 

0  .  L  2  0 

ENTER  THE  SECOND 

ROW. 

A  ■  A  A  /  \ 

!  V  •  A  /W\ 

'/  V  '  •  v 

1  .  9  00 

0 . 0  (  ',  0 

9 .  nor. 

ENTER  T 

HE  THIRD 

ROW. 

X  .  X  X  X 

X .  XXX 

A  .  A  A  A 

0 . 000 

-0.L2S 

0.906 

TRANSFORMS T ION  MATRIX 


C!  C2  CO.  IN  THIS  FORMAT: 

CL  CS  CS.  IN  THIS  FORMAT 
C7  C3  CO.  IN  THIS  FORMAT: 
FOR  DIPOLE*  3 


CL  CS  C6 
C7  CS  C 3 

ENTER  THE  FIRST  ROW.  C!  C2  C3.  IN  THIS  FORMAT: 

X.XXX  X.XXX  X.XXX 

t . non  0.9  o  n  n . n o o 

ENTER  THE  SECOND  ROW.  CL  CS  CG.  IN  THIS  FORMAT 

\/  \/  V  »/  ./V«/  *./  «/  '/ 

A  •  <  v  A  A  /\  .  mhA  /•  >  i\ /\ /' 

0.9  no  l.OC.n  9.009 

ENTER  THE  THIRD  ROW.  C7  C3  C9.  IN  THIS  FORMAT: 
X.XXX  X.XXX  X.XXX 
0.000  0.000  1.000 

ENTER  THE  MEAS  RADIUS  (1ST  GRAD  COIL):  X.X  <M) 


Appendix  B 


This  section  contains  the  IBM  30/30  control  language 
used  to  run  the  Surface  II  Graphics  plot  program.  The 
software  consists  of  a  variety  of  routines.  A  desired  plot 
can  be  formed  by  listing  the  pertinent  routines  and 
specifying  their  reguired  parameters  in  a  parameter  file.  A 
detailed  explanation  of  the  capabilities  of  each  routine  and 
how  to  specify  their  parameters  is  given  in  reference  13. 


for  XYJPLQI 


f.i f r  • i  '  /hi.  V  •['  '■  l  •  T  T  .  dE.VE  .  tt'M-.'Li.iT  _  -  A  .  ,'i<j  r  .V  <  M :  fvO .  1 1, IE-  (  3) .  CLASS-F 
. i o ( i l  : [:  dm  /  p.-.. :  . surf ;• . .. :>pdl i e* . r* i hr 

. ;  1  i  f.l  a :  i  1 ..  T  0  •;  i.1  •  [>Lt  I  -  :  '  MP  *  .  UNI  t>h  1  :  I ' .  SPACE  •  'L  •  .  •  ?.!.! .  :.0  )  ) 

,101.1:1'  .  r  ■  ii^jr  :!!)  •  L)  I>  Licit  :  iE.'-IPL  .  ONI’  =  0 1  '-7 .  SPACE  *  i  i.M-  .  ;  5n  .  3U  )  ) 

,i0 OF ;  1  //:  ’oOFo-ii  dlj  Dc.;i  =  t  iv: .  on :  i .  [>  i  i.v .  spas F.  »  <  TR.-: . <  Su . 30  > » 

.10070  /  / F T 0 :  F 0 0  !  DD  D M- ~F  TP  .■[> .  :  It  -  YPL 0  T )  .  Di 3P-  3.-.R 

,1000!!  ,  /-rooroo  l  DD  :'.YT.0:r  -A 

00  030  //F  TO/FOU  r  D 0  6YS3UT-A 

00  1  00  //FT  T  IF  00  i  DD  DSN-HE7PJD.  oYF'^OF  .[)  I  :  p  -  5  H  R 

00110  //FT  1  OF 00  1  DD  DSN-HF TP JD.PL01 ( OUT.  INE ) . OISP-SHR 

00  1  20  //FT  i  F  F  0  0  1  DD  DSUm.vSRFGUT  1  ..JNIT-DiSK. SPACE  m  TRK.  (2.2)  ) 

00130  //FT 2 OF  00  1  DD  DSN*  <t.*SRF3UT  2 . IJN I T*D ISK . SPACE*  <  TRK .  (2.2)  ) 

00  100  / / F T o 2 F 0 0 1  DD  SYSOUT-A 

00 150  //FTO3F001  DD  SYSGUT-A 

00160  //FT3IF00!  DD  SYSOUT-A 

00170  / / * F T 3 1 F 0 0 1  DD  DSN- JC0RJGJ .SRF0UT3 . DATA .UNIT-DISK. SPACE- ( TRK, ( 1 , 1 )) , 
001  SO  //•*  D I  j P  =  f  .KEEP)  . DCB- ( BLKSI ZE-800 , LRECL-800  . RECr M=U ) 

00190  //PLOTLOG  DD  SYS0UT-A 

00200  //VECTRI  DD  DSN-HEGBP . DSG. VECTR1 .UNI T-DISK , SPACE- ( TRK. (2 .2) ) , 
00210  //  DISP- ( OLD . KEEP ) 

00220  //VECTR2  DD  DSN-HEGBP . DSG . VECTR2  .UNIT-DISK . SPACE- ( TRK. < 2.2 >) , 
00230  //  DI3P- ( OLD . KEEP ) 

0024D  //PL3TPARM  DD  DSM-HETPJD • P10TPARI1 . DATA ( PLOTDSG 1 ) . DISP-SHR 

00250  // 

END  OF  DATA 
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AD-A152  227 
UNCLASSIFIED 


A  NULTIPOLE  MODEL  OF  THE  OBSERVED  CEREBRAL  CORTEX 
HAGNETIC  FIELD(U)  AIR  FORCE  INST  OF  TECH 
MRIGHT-PATTER50N  AFB  OH  SCHOOL  OF  ENGINEERING 
P  J  DE  REGO  DEC  84  RFI T /GE/ENG/84D-23  F/G  6216 


MICROCOPY  RESOLUTION  TEST  CHART 

fs| A 1  ><  >N A|  HIlWA"  ■*■ 


» 


I 


I 


I 


00010 
0  G  0  2  0 
00030 
00040 
000S0 
00060 
00070 
0  C  0  8  0 
00100 
00110 
00120 
or,  1 30 
00140 
00150 
00  1-GO 
END  OF 


TITLE 

DEVICE 

IDXY 

ROUTLINE 

ECHO 

GRID 

EXTREMES 

BLANK 

CONTOUR 

CINTERVAL 

POUTLINE 

SIZCONTGUR 

BOX 

PERFORM 

STOP 

DATA 


MAGCAL  -  CONTOUR  PLOT  USING  CALCULATED  DATA-SPHROS 
6.  'LT  DE  REGO' 

170. 1 1 .3.1 .2.3.-1 .0. I ,0, ’ (2F3.3.SX.EI5.8) ' 

12. 1 . ' ( 10X.2F10.2) ' 

0. SO. SO. 0.0. 0.0 

-12.12,-12.12 


0.0.0.2.0.S.0, 1 .3.0. 1 .10 

1.6.6 

1.4. 1.4. 0.0, 0.1, 0.1 


far.  PLat 


00010  //i-IETPJD;.’  JOB 
00015  //  CLASS'D 


00020  //PLOT  c 
00030  //STEFLIB 
00040  //PLOTLOG 
OGOSO  //VECTR1 
00060  //VECTR2 
00070  //SYSVECTR 
00080  // 


■1  HE T .  DEREGO .  NO TIFY*HET P JD .  T I.'IE  •  ( 2 ) ,  MSGCLASS  *A . 

1C  PGM =VT PLOT , COND* ( 0 .LT ) 

DD  DSN=VTEC.PLOTLIB.DISP*SHR 
DD  SYSOUT'A 

DD  DSN'HEGBP . DSC- .  VECTR1  . D ISP »( OLD . KEEP  > 

DD  DSN=HEGBP . DSG . VECTR2 . D ISP- ( OLD . KEEP > 

DD  UN  I T  -  CEE . DCS  = ( BLKS I ZE  =  1 33 . LRECL  =  133. RECFM- F  A ) 


Appendix  C 


This  section  contains  the  data  generated  by  the  program 
MAGCAL2  for  each  of  the  five  models.  The  same  dipole 
parameters  were  used  in  the  actual  experiments.  The  data 
collected  from  the  experiments  is  given  in  appendix  D. 
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0 . 443972SSE+00 

BD  ( 

3,15.0.1098) 

-6.674 

-1  . 323 

8.584 

0.333/42S4E+0G 

BD  ( 

3.16.0.1 093 ) 

-6 . 168 

-2.970 

8.534 

0.21847S76E+00 

BD  ( 

3.1 7.0.1098) 

-5.352 

-  4 . 263 

6.584 

0 . 66254  7  76E -0 1 

BD  C 

3.18.0.1 098 ) 

-4.268 

-5.352 

8 .534 

-0. 50292 108E-01 

BD  ( 

3.19.0.1 098 ) 

-2.970 

-6. 168 

6.534 

-0. 18431 637E+00 

BD( 

3.20.0. 1093) 

- 1 .523 

-6.674 

8.534 

-0 . 30909926E  +00 

BD  ( 

3.21 .0.1098) 

C.OOC 

-6.846 

6.534 

-0 . 41 836264E+00 

BD  f 

3.22,0.1098) 

1.523 

-6.674 

8.534 

-0 .506636 1 SE  +  00 

BD  ( 

3,23.0.1098' 

2.97C 

-6.166 

6.534 

-  0. 569562 16E +  00 

BD  C 

3.24.0. 1098) 

4.268 

-5.352 

8.584 

-0 . 6039 1 653E  +00 

BD  f 

3,25.0.1093) 

5.352 

-4.266 

6.534 

-  0 . 60  796809E  +  00 

BD  ( 

3.26.0. 1098) 

6.168 

-2.970 

3.534 

-0 . 58 1 S3307E+00 

BD  ( 

3,27,0.1098) 

6.674 

-1 .523 

8.534 

-0 .52593303E+00 

BD  ( 

4.  0,0.1098! 

8 . 535 

0  .  000 

6 . 846 

-0.31 1 35607E+00 

BD( 

4.  1.0.1093) 

6.369 

1  . 9 1  C 

6.846 

-0 . 25651 99  7E  +00 

BD( 

4,  2,0.1098) 

7.734 

3.725 

6.846 

-0. 18882084E+00 

BD  ( 

4,  3.0.1098) 

6.712 

5.352 

6.346 

-0.1 1 1 65357E+00 

BD  ( 

4,  4.0. 1098) 

5.352 

6.712 

6.846 

-0 . 2888755 1 E-0 1 

BD  ( 

4.  5.0.1098) 

3.725 

7  .  734 

6.846 

0.35327006E-01 

BD( 

4.  6.0.1U98) 

1  .910 

8.369 

6.346 

0. 13676739E+00 

BD< 

4.  7.0.1098) 

-C.OOC 

8.535 

6.846 

0.21 134979E+00 

BD( 

4,  8,0.1098) 

-1.910 

8.369 

6.846 

0 . 27533382E+00 

8D( 

4.  9,0.1093) 

-3.725 

7.734 

6.846 

0.32551 146E+00 

I 
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Calculated  Eaia  Hotel  M 


DIPjL:»! 

DIF’GLE>2 

DIF0LE*j 

MLMEM 

i  11%%' -0 

c 

0 . 

55E-05 

0.5SE-U5 

i.jCAT: 

1  o.Obb 

0. 

053 

0.040 

i  rvrl!  J  :■ 

: .  ■  i  ii  o  o 

.00  0  0 . 

0 o (I  0  . 

000  1.000 

0.000  0.000  1.000 

MATRi; 

H  11.000  1 

.000  0. 

non  - 1 . 

0 0 Li  0.000 

0.000  -1.000  0.000 

0.000  0 

.000  1. 

000  0. 

000  0.000 

1  .000  0.000  0.000 

THE  'IE AS  SF'HERL  RADIUS  IS 

0.  1098 

M 

Sr  HER 

DATH  POINTS : 

THETA=T  *0.2244: 

PHI =F  *0 . 

2244 

BD(T.F.R) 

A 

Y 

Z 

FLUX  DENS.  XEE-11 

BD  ( 

0.  0.0. 1 098 ) 

0 . 0  0  0 

n.000 

•  10.980 

O.OOOOOOOOE+OO 

BD  < 

1  .  G.O. 1093) 

2 . 445 

0.000 

•  10.705 

-0 . 424 1  0666E  +  00 

BD( 

1.  1.0.1098) 

o  '•  o  O 
C.  •  JUL 

0.544 

10.705 

-0.2141 3589E+00 

BD  C 

1  .  2.0.1098) 

2.201 

1 . 060 

10.705 

0 . 6S727420E-02 

BD  f 

1,  3.0.1098) 

1.910 

1 .523 

10.705 

0 . 22695 1 66E+00 

BD  C 

1  .  A. 0. 1093) 

1  .523 

1 .910 

10.705 

0 . 43594374E+00 

E'D  ( 

1.  S. 0.1093) 

1  .060 

2.201 

10.705 

0 . 62308735E  +  00 

BD  ( 

1.  S, 0.1093) 

0.544 

2.382 

10.705 

0 . 778981 33E+00 

BD  C 

1.  7.0.1098) 

-C.00C 

2 . 443 

1C.  705 

0 . 89581 430E+00 

BD  ( 

1  ,  8.0.1098) 

-0.544 

2 . 332 

10.705 

0 . 967  725S3E  +  00 

BD< 

1  .  9.0.1098) 

-1 .060 

2.201 

10.705 

0.991 1 1 1S3E+00 

BD  ( 

1  .10.0.1098) 

-1.523 

1.910 

10.705 

0 . 96479923E+00 

BD  ( 

1.11.0.1098) 

-1 .910 

1  .523 

10.705 

0.8901 C305E+00 

BD! 

1  .  12.0. 1098) 

-2.201 

1  .060 

10.705 

0.770 7817SE +00 

BD  ( 

1.13.0.1 093 ) 

-2.382 

0.544 

1C. 705 

0.61 28C584E+00 

BD  f 

1.14.0. 1098) 

-2.443 

-0.000 

10.705 

0.42410010E+00 

BD  f 

1  .15.0. 1098) 

-2.332 

-0.544 

10.705 

0.21 41 2909E+00 

BD  f 

1 . 18,0.1 098  ) 

-2.20' 

-1.060 

10.705 

-0.65801 106E-02 

BD  ( 

1  .1  7.0.1  098) 

-1  .910 

-1  .523 

10.705 

-0 . 22695762E +00 

BD  C 

1.18.0.1 098 ) 

-1  .523 

-1.910 

10.705 

-0 . 4359S370E+00 

BD  ( 

1.19.0.1 093 ) 

-1 .060 

-2.201 

10.705 

-0.62309307E+00 

BD  C 

1.20.0. 1098) 

-0.544 

-)  OA') 
J  .  •  j  C. 

10.705 

-0 . 77898550E+00 

BD  C 

1 .21 .0.1093) 

0.000 

-2.445 

10.705 

-0.89531 704E+00 

BD< 

1.22.0. 1098) 

0.544 

-2.332 

10.705 

-  0 . 967  7273SE  +00 

BD' 

1 . 25 .0.1 093 ) 

1  . 060 

-2.201 

10.705 

-0.9911 1 187E+00 

BD  < 

1,24.  IT  .  I  098) 

1  .523 

-1.91  0 

10.705 

-  0 . 964  79762E  +00 

BD' 

1  .25.0. 1098) 

1.91  0 

-1  .525 

10.  705 

-0 . 690 1 G495E+00 

BD' 

t . 2b .0.1 098 ) 

2. 20  1 

-  I  .  060 

10.705 

-0.7707771  IE  +  00 

BD  ( 

I  .27.0. 1  093) 

v>  " 

-i,  .544 

11..  7D5 

-0 . 6 1  28 000 BE  +  0U 

E:[)< 

2.  0  .  U . ' 0 8 8 ) 

4  .  /'  6  4 

|| .  000 

9 . 89  : 

-n .  59.89/.:y  7E-h.hi 

BD  ( 

2.  1.0.1093) 

4.645 

1  . 060 

S .  38 3 

-0  .  35611500'!  E  +  1II.1 

BD! 

2.  2,0.1098) 

4.292 

2.067 

9 .  ;.i3  j 

-0. 1 55o46A2E ‘DU 

BD  f 

2.  5.0.1093) 

3.725 

2 .9  70 

2-  •  O  3  ^ 

0.5  TT8'.*46Rr  Ii  ! 

BD! 

2.  4.0.1098) 

2.97  0 

7’C 

9.393 

n .  25905 63 3 E  *0o 

BD! 

2.  5.0.1093) 

<1.06/ 

4 . 292 

9.893 

0 . 4bi  ''3256I  *-i.iii 

BD! 

2.  6.0.1093) 

1  . 060 

4.645 

9.893 

0 . 622 1 4845E *0o 

BD  ( 

2.  7.0.1093) 

-o.oor. 

4.  764 

3.393 

0.7611 668 7E *00 

BD! 

2.  8.0.1098) 

-1 .060 

4.645 

9.893 

0.362U 1 6 1 % *00 

BD! 

2.  9.0.1093) 

-2.067 

4.232 

5.893 

0.91  9640241.  *00 

BD! 

2.10.0.1098) 

-2.970 

3.725 

9.893 

0.931 1 4996E  +  00 

BD! 

2.11.0.1 093 ) 

-3.725 

2.370 

9.333 

0.895%  7  78E +00 

BD! 

2.12.0.1098) 

-4.292 

2.067 

9.893 

0.81 535747E  +00 

BD! 

2.13.0.1098) 

-4.645 

1  .060 

5 . 333 

0.69483703E+00 

BD! 

2.14,0.1098) 

-4.764 

-0.000 

9 . 893 

0 . 5389732 1 E  +  00 

1%,'P+M) 

Cl  r  R  Orv'N  ) 

0.000 
0.000 
1  .000 


TESLA 


BD( 13,22,0. 1098)  0.944  -2.382  -10.709  -0  . 92279234E- 02 

BD (  1 3 . 23. 0.1098)  1.0S0  -2.201  -1  0.709  -0 . 994420 1 0E - 02 

BD( 13,24,0. 1098)  1.923  -1.910-10.709  -0 . G3827998E- 02 

80(13,29.0.1098)  1.910  -1 .923  -1 C . 709  -0 . 64621 344E-02 

BD< 13,26,0. 1098)  2.201  -1.060-10.709  -0 . 62374324E-02 

60(13,27,0.1098)  2.382  -0.944  -10.709  -0 . 96999949E-02 

BD < 14.  0,0.1098)  -0.000  0.000  -10.980  - 0 . 69667046E- 08 

END  OF  DATA 
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BDM  1  ,  12.0. 1038) 

-6. 168 

2.970 

-8.585 

0.21 292880E-0 1 

6D( 1 1 ,13.0.1008) 

■6.674 

1  .523 

-3.585 

0. 19457854E-01 

BD ( 11.14,0. 1098) 

-6.846 

-0.000 

-8.585 

0.166471 15E-01 

BD( 1 1 . 15.0. 1098) 

-6.674 

- 1 . 523 

-8.585 

0.130C1624E-01 

BD  < 1 1,18,0. 1098) 

-6. 168 

-2.970 

-8.585 

0 .8704 1 594E-02 

80(11 .17,0.1098) 

-5.352 

-4.268 

-8.585 

0.397C2505E-02 

80(11,18,0. 1098) 

-4.263 

-5.352 

-8.585 

-0.96275867E-03 

80(11 ,19,0.1098) 

-2.970 

-6.168 

-8 . 565 

-0 . 5847491 3E-02 

BD ( 1 1 ,20,0.1098) 

-1.523 

-6.674 

-8.585 

-0. 10433986E-01 

BD ( 1 1 ,21 ,0. 1096) 

0.000 

-6.846 

-8.585 

-0. 145C7055E-01 

BO ( 1 1 ,22,0. 1098) 

1.523 

-6.674 

-8.585 

-0. 17847657E-01 

BDM  1  ,23,0.1098) 

2.970 

-6.168 

-8.585 

-0.20293307E-01 

BD ( 11,24.0.1098) 

4.268 

-5.352 

-8.585 

-0.21721 352E-0 1 

80(11,25,0.1098) 

5.352 

-4.268 

-8.585 

-0 . 220601 93E-0 1 

BD ( 1 1 .26.0. 1098) 

6.  168 

-2.970 

-8.585 

-0 . 21 292843E- 0 1 

BD ( 1  1  ,27.0. 1  098) 

6.674 

-1.523 

-8.585 

-0. 19457784E-01 

BD ( 12,  0,0.1098) 

4.764 

0.000 

-9.893 

-0. 1 02324 1 5t -0  1 

BD ( 1 2 .  1.0.1098) 

4.645 

1 .060 

-9.893 

-0.79920925E-02 

BO ( 12,  2,0.1098) 

4.292 

2.067 

-9.893 

-0 . 535 1 0033E  -02 

BD ( 12,  3,0.1098) 

3.725 

2.970 

-9.893 

-0 . 2441 591  4E  -  02 

BD ( 12,  4,0.1098) 

2.970 

3.725 

-9.893 

0.59024827E-03 

BD ( 1 2 .  5,0.1098) 

2.067 

4.292 

-9.893 

0.35924902E-02 

BD ( 12,  6,0.1098) 

1.060 

4.645 

-9.893 

0.641 4588SE-02 

BD ( 1 2 ,  7,0.1098) 

-o.ooc 

4.764 

-9.893 

0.691 SG406E  -02 

BD ( 12.  8,0.1098) 

-1.060 

4.645 

-9.893 

0.  10968439E-01 

BD ( 1 2 ,  9,0.1098) 

-2.067 

4.292 

-S.893 

0.12471 83BE-01 

BD( 12, 10,0. 1098) 

-2.970 

3.725 

-9.893 

0. 13349853E-U1 

BD ( 1 2 , 1 1 ,0.1098) 

-3.725 

2.970 

-5.893 

0 . 1 3556436E  -0 1 

BD( 12. 12,0. 1098) 

-4.292 

2.067 

-9.893 

0. 13087146E-01 

BD ( 1 2 , 13. 0.1098) 

-4.645 

1.060 

-9.893 

0. 1 1959597E-01 

BDM2, 14,0. 1098) 

-4.764 

-0.000 

-9.893 

0. 10232352E-01 

BD< 12, 15,0. 1098) 

-4.645 

-1 .06C 

-9.893 

0 . 799201 80E-02 

BD ( 12, 16,0. 1098) 

-4.292 

-2.067 

-9.893 

0.5350S139E-02 

80(12,17,0.1098) 

-3.725 

-2.97G 

-9.893 

0 . 2441 50  78E  - 02 

BD ( 12.18.0.1098) 

-2.970 

-3.725 

-9.893 

-0 .5903365 1 E -  03 

BD(1 2, 19, 0.1 098) 

-2.067' 

-4.292 

-9.893 

-0 . 3592581 3E  -02 

BD ( 12.20,0. 1098) 

-1.060 

-4.645 

-9.893 

-0.641 46742E-02 

BD ( 12,21 ,0.1098) 

C.OOC 

-4.764 

-9.893 

-0.69151 077E-0 2 

BD ( 12.22,0. 1098) 

1  .060 

-4.645 

-9.893 

-0.  10968495E-01 

BD ( 12.23,0. 1098) 

2.067 

-4.292 

-9.893 

-0 . 1  24  7 1  88 1  E -  0 1 

BD ( 12.24,0. 1098) 

2.970 

-3.725 

-9.893 

-0. 13349872E-01 

BD( 12.25,0. 1098) 

3.725 

-2.970 

-9.893 

-0. 13556425E-01 

BD ( 12.26.0. 1098) 

4.292 

-2.067 

-9.893 

-0. 130871 16E-01 

BD(12. 27, 0.1098) 

4.645 

-1  .060 

-9.893 

-0. 1 1955553E-01 

B D ( 13,  0,0.1098) 

2.443 

0.000 

-10.705 

-0.4876721 7E-02 

BD ( 13.  1.0.1098) 

2.382 

0.544 

-10.705 

-0.38069016E-02 

BDM  3.  2,0.1098) 

2.201 

1  .060 

-10.705 

-0 . 25500872E -  02 

BD ( 1 3 .  3.0.1098) 

1  .310 

1  .523 

-10.705 

-0.1 1633986E-02 

B  D  < 13,  4,0.1098) 

1  .523 

1.910 

-10.705 

0.281 623 1 0E-03 

BD ( 1 3  .  5,0.1096) 

1 .060 

2. 201 

-10.705 

0.1 7 1 25262E  -02 

BD (  13,  6,0.1098) 

0.544 

2.382 

-10.705 

0 . 30S75558E- 02 

BD ( 1 3 ,  7,0.1098) 

-0.000 

2.443 

-10.705 

0.42492673E-02 

B D ( 13,  8,0.1098) 

-0.544 

2.382 

-10.705 

0.52278973E-02 

BD( 1 3 ,  9,0.1098) 

-1.060 

2.2C1 

-10.705 

0.594433246-02 

80(13,10,0. 1098) 

-1.523 

1.910 

-10.705 

0.63627884E-02 

BDM3.1  1  ,0.1098) 

-1  .910 

1  .523 

-10.705 

0 . 64621 381 E -02 

BD( 13. 12,0.1098) 

-2.201 

1  .060 

-10.705 

0.62374435E-02 

BD(13, 13. 0.1098) 

-2.382 

0.544 

-10.705 

0 . 569997  73E-02 

BDM  3,  14,0.  1098) 

-2.443 

-0.000 

-10.705 

0 . 487669 1 8E- 02 

BDM3. 15.0.1098) 

-2.382 

-0.544 

-10.705 

0.38088665E-02 

BDM  3,  16,0.  1098) 

-2.201 

-1.060 

-10.705 

0.25500439E-02 

BD< 13, 17,0. 1096) 

-1 .910 

-1  .523 

-10.705 

0.11  633583E -  02 

BDM3. 18,0.  1098) 

-1.523 

-1.910 

-10.705 

-0.281 66780E -03 

BDM3, 19.0. 1098) 

-1 .060 

-2.201 

-10.705 

-0. 1 7125721E-02 

BDM3. 20,0.  1098) 

-0.544 

-2.382 

-10.705 

-0.30575949E-02 

BD( 13,21 ,0.1098) 

0.000 

-2.443 

-10.705 

-0.42493008E-02 
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BD<  9,  2.0.1098) 
BD (  9,  2.0.1098) 
BO (  9,  4,0.1 098 ) 
BD C  9.  5.0.1098) 
BD(  9,  6.0.1098) 
BD f  9.  7.0.1098) 
BD !  3.  8,0.1098) 
BD C  9.  9,0.1093) 
ElD <  3,10,0.1098) 
BD f  9.1  1.0.1098) 
BD (  9,12,0.1098) 
BD (  9.12.0.1098) 
BD  C  3,14,0.1098) 
BD <  9.15.0.1098) 
R[)(  9.16.0.1098) 
BD(  9.17.0.1098) 
BD (  9.18.0.1098) 
80 (  9,19,0.1098) 
BD C  9.20,0.1098) 
BD (  9,21,0.1093) 
BD C  9.22,0.1098) 
BD C  9.23,0.1098) 

BD <  9,24,0.1098) 
BD (  9,25.0.1098) 
BD(  9,26,0.1098) 
BD<  9.27,0.1098) 
BD (10,  0,0.1098) 
BD ( 1 0 ,  1,0.1098) 
BD ( 10,  2,0.1098) 
BD ( 10,  3,0.1098) 
BD (10,  4,0. 1098) 
BDMO,  5,0.1098) 
B D (  10,  6,0.1098) 
BDMO,  7,0.1096) 
BD ( 1 0 ,  8,0.1098) 
BD ( 1 0 ,  9,0.1038) 
BDMO,  10,0.  1038) 
BD  MO,  1  1  ,0.  1098) 
BD ( 10, 12,0. 1098) 
BD ( 10, 13,0. 1096) 
BD( 10,14,0. 1098) 
BD110, 15,0. 1098) 
BD<  10, 16,0. 1098) 
BDMO, 1  7, 0.1098) 

BDMO,  18,0.  1098) 
BDMO, 19,0. 1098) 
BDMO, 20,0.  1098) 
BDMO, 21 ,0.1098) 
BDMO,  22,0.  1098) 
BDMO,  23,0. 1098) 
BD ( 10,24,0.1098) 
BDMO,  25, 0.1  098) 
BDM  0,26,0.  1098) 
BDMO, 27, 0.1  098) 
BDM  I  ,  0,0.1098) 
BDM  1  ,  1,0.1  098) 
BDM1  ,  2,0.1098) 
BD ( 1 1 ,  3,0.1098) 
BDM  1  ,  4,0.  1098) 
BDM  1  ,  5,0.1098) 
BDM  1  ,  6.0.  1098) 
BDM  1  ,  7.0.1093) 
BDM  1  .  8,0.1098) 
BDM  1  ,  9,0.  1098) 
BDM  1  ,  10,0.  1098) 
BDM  1  .11  .0.1  098) 


8.913 

4.292 

-4.764 

-0. 18944908E-01 

7.734 

6 . 1 63 

-4.764 

-0.858010726-02 

6.168 

7.734 

-4.764 

0.221493966-02 

4.292 

S.913 

-4.764 

0.12898933E-01 

2.201 

9.645 

-4.764 

0.229361 176-01 

C.OOC 

9.893 

-4.764 

0.318232106-01 

2.201 

9.645 

-4.764 

0.391 145016-01 

4.292 

8.913 

-4.764 

0.444444496-01 

6.168 

7.  734 

-4.764 
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BDf  3.12,0.1098) 
BDf  3,13,0.1098) 
BDf  3,14.0.1098) 
BDf  3,15.0.1098) 

BDf  3.16.0.1098) 
BDf  3.17,0.1036) 
BDf  3.18,0.1093) 
BDf  3.19,0.1096) 
BDf  3.20,0.1098) 
BDf  3.21 ,0.1096) 
BDf  3,22.0.1098) 
BDf  3.23,0.1096) 
BDf  3,24,0.1093) 
BDf  3.25,0.1098) 
BDf  3,26,0.1098) 
BDf  3,27,0.1098) 
BDf  4,  0,0.1098) 
BDf  4,  1,0.1096) 
BDf  4,  2,0.1098) 
BDf  4,  3,0.1096) 
BDf  4,  4, U. 1098) 
BDf  4.  5,0.1096) 
BDf  4,  6,0.1098) 
BDf  4,  7,0.1096) 
BDf  4,  8,0.1098) 
BDf  4,  9,0.1098) 


-4 
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-1  . 

060 

9. 

893 

-4 

.292 

-2. 
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9. 
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-3 
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9. 

893 

-2 
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-3. 
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9. 
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-2 

.067 

-4. 
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9. 

893 

-1 

.060 

-4. 

645 

9. 

893 

0 

.000 

-4. 

764 

9. 

893 

1 

.060 

-4. 

645 

9. 

893 

2 

.067 

-4. 

292 

9. 

893 

2 

.970 

-3 . 

725 

9. 

893 

O 

.725 

9  70 

9. 

893 

4 

.292 

-2. 

067 

9. 

893 

4 

.645 

-1  . 

060 

9. 

893 

6 

.846 

0. 

000 

8. 

584 

6 

.674 

1  . 

523 

8. 

534 

6 

.  1 68 

2. 

970 

8. 

584 

c 

.352 

4. 

268 

8. 

584 

4 

.268 

5. 

352 

8. 

534 

2 

.970 

6. 

168 

8. 

584 

1 

.523 

6. 

674 

8. 

584 

-C 

.000 

6. 

846 

8. 

584 

-1 

.523 

6. 

674 

8. 

584 

-2 

.970 

6. 

1  68 

8. 

584 

-4 

.268 

5. 

352 

8. 

584 

-5 

,352 

4. 

268 

8. 

534 

-6 

,  168 

2. 

970 

8. 

584 

-6 

,674 

1  . 

523 

8. 

584 

-6 

.846 

-0. 

000 

8. 

584 

-6 

.674 

-1  . 

523 

8. 

584 

-6 

.  163 

-2. 

970 

8. 

584 

-5 

.352 

-4  . 

268 

8. 

564 

-4 

.268 

-5. 

352 

8. 

584 

-2 

.970 

-6. 

168 

8. 

584 

-T 

.523 

-6. 

674 

8. 

584 

0 

.000 

-6. 

846 

8. 

564 

i 

.523 

-6. 

674 

8. 

584 

2 

.970 

-6. 

168 

8. 

564 

4 

.268 

-5. 

352 

8. 

584 

5 

.352 

-4. 

263 

8. 

5S4 

6 

.  168 

-2. 

970 

8. 

584 

6 

.674 

-1  . 

523 

8. 

564 

8 

.585 

0. 

000 

6. 

846 

8 

.369 

1  . 

910 

6. 

846 

7 

.  734 

3. 

725 

6. 

84b 

6 

.712 

5. 

352 

6. 

846 

c 

U 

.362 

6. 

712 

6. 

846 

3 

,  725 

7. 

734 

6. 

846 

1 

,910 

8. 

369 

6. 

846 

-0 

,000 

8. 

585 

6. 

846 

-1 

.910 

8. 

369 

6. 

846 

-3 

,725 

7. 

734 

6. 

846 

0 . 2831 6540E  +00 
0.61 380204E-0 1 
-0 . 1 63431 S9E+0C 
-0 . 380 1 471 OE  +  OO 
-0 . 57  774985E  +  00 
-0 . 746381 94E+00 
-0 .87758696E+0C 

-0 . 96478570E+00 
-0. 100360S8E+01 
-0 . 992 1 0 1 S5E+00 
-0 . 930848 1 8E+00 
-0 . 82291 776E+00 
-0 . 67372340E+00 
-0 . 397 1 891 4E+00 
-0.26899797E+00 
-0. 12731850E+00 
0 . 20  74SS08E-0 1 
0. 16776884E+00 
0 .30637980E+00 
0 . 42962736E+00 
0.531 3321 4E  +  00 
0 . 606393 1 OE+OO 
0.651 04681 E  +  00 
0 . 66305465E+00 
0 . 64 1 81 352E+00 
0 . 58838892E+00 
0 .50S46032E+00 
0 . 397 1 8533E+00 
0 . 26899403E+00 

0. 1 2731 36 7 E  +  00 
-0.2074971  IE-01 
-0. 16777325E+00 
-0 . 30638385E+0C 
-0 . 42963099E+00 
-0.53133494E+0C 
-0 . 606395 1 8E+00 
-0.651 0480CE+0G 
-0 . 66305435E+00 
-0.641 81 232E+0G 
-0 . 5883867 1 E+00 
-0.S0S4S728E+0G 
-0.27379555E+00 
-0.1 997231 2E+00 
-0. 1 1563587E+00 
-0.2S750104E-01 
0.65426707E-01 
0.15332288E+0G 
0.23353070E+00 
0 . 3020286GE+00 
0.355381 19E+00 
0.39091331E+0G 
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Data  for  Model  11 


DIPOLE'/ 1 

MOMENT  0. SEE -OS 

LOCO  IN  il .  066 


D IP0Lc#2 
0.55E-05 
0 . 053 


DIP0LE#3 

0.55E-05 

0.040 


trams 

MATRIX 


t  .  IJ  0  0  0 . 0 1.1 0  0 .000 

o.ooo  i .  o  i.  n  o.n  n  o 
o.ooo  o.ooo  l.orm 


0.000  1.000  0 
-1.0011  0.000  0 
o.ijoo  o.ooo  i 


000  0.707  0.707 

OHO  -0.70/  0.707 

000  0.000  u.ooo 


IHE  ME AS  SPHERE  RADIUS  IS  0.1098  M 


a  *• 

«a 


HER 

DA  I A  POINTS: 

THE  TA  =  T  *0 

.2244; 

PHi=r»u 

2244 

HD ( T . F . R ) 

\/ 

/\ 

Y 

z 

FLUX  DENS.  XEE- 

Bin 

0.  0.0.1098) 

0 .  ooo 

0.000 

10.980 

O.OOOOOOOOE+OO 

BDl 

1.  0.0.1098) 

2.443 

0.000 

10.705 

-0.391 SC470E+00 

BD 

i.  1.0.1098) 

2 . 382 

0.544 

10.705 

”0. 1 6483694E+00 

BDl 

1.  2.0.1098) 

2.201 

1 .060 

10.705 

0.70096433E-01 

BD  C 

1,  3.0.1098) 

I  .  9 1  0 

!  .  523 

10.705 

0 . 30 ! 51 480E+00 

BDl 

1  .  4,0.1 098 ) 

1  . 523 

1  .910 

i  c :  7  u  s 

0.51 781 327E+00 

BD  ( 

1  .  5.0.1098) 

1  .060 

2.201 

10.705 

0.7081 4663E+00 

BDl 

1  .  8,0.1098) 

0.544 

2.382 

10.705 

0 . 86297077E+00 

BDl 

1.  7.0.1098) 

-c.ooc 

2.443 

10. 705 

0 . 97452235E+00 

BDl 

1.  8,0.1098) 

-0.544 

A_  .  JO  4. 

10.705 

0. 1 0372057E+01 

BDl 

1  .  9.0.1 098 ) 

-1 .060 

2.201 

10.705 

0 . 1 04  78732E+0 1 

BDl 

1,10.0. 1098) 

-  1 . 523 

1.910 

10.705 

0 .  J  0060081 E+0 1 

BDi 

1.11,0.1098) 

-1 .910 

1  ,  S23 

1C.  705 

0.91 36921 8E+00 

BDl 

1,12.0.1098) 

-2.201 

1  .060 

10.705 

0.7755581 1 E+00 

BDl 

1  .13.0.1098) 

-2.382 

0.544 

10.705 

0.59853494E+00 

BD( 

1,14,0.1 098 ) 

-2.443 

-0.000 

10.705 

0.391 49/5SE+00 

BDl 

1 ,15,0.1098) 

-2.382 

-0.544 

10.705 

0. 16482967E+00 

BDl 

!.  16.0. 1098) 

-2.201 

-1.060 

10.705 

-0.701041 22E-0 1 

BDl 

1,17,0.1098) 

-1.910 

-1 .523 

10.705 

-0.301521 06E+00 

BDl 

1,18,0.1098) 

-1 ,523 

-1.910 

10.705 

-0.51781 940E+00 

BDl 

1  ,19,0.1098) 

-1  .060 

-2.201 

10.705 

-0.7  081 5235E  +  U0 

BDl 

1.20,0. 1098) 

-0.544 

-2,382 

10.705 

-0 . 86297494E+00 

BDl 

1  ,21 ,0.1  098) 

O.OOO 

-2.443 

10,705 

-0.97452497E+0U 

BDl 

1,22,0. 1098) 

0.544 

-2.382 

10.705 

-0. 10372066E+01 

BDl 

1  ,23,0.1098) 

1  . 1)60 

-2.201 

10.705 

-0. 10476783E+01 

BDl 

1,24,0.1098) 

1  .523 

-1,910 

10.705 

-0. 10060062E+01 

BDl 

1  ,25,0.1098) 

1.910 

-1 .523 

1C.  705 

-0 . 91 3688S4E+00 

BDl 

1,26,0. 1098) 

2.201 

-1.060 

10.705 

-0 . 77555287E+00 

BDl 

1,27,0. 1090) 

2 . 382 

-0.544 

10.705 

-0.59852856E +00 

BDl 

2.  0,0.1098) 

4.764 

0.000 

9.893 

-0 . 49075735E+00 

BDl 

2,  1,0.109b) 

4 .645 

1  .060 

9 . 853 

-0 . 283 1 720/E+0G 

BDl 

2,  2,0.1098) 

4 . 292 

2.067 

9.893 

-0.61 387379E-0 1 

BDl 

2,  3,0.1098) 

3.725 

2.970 

9.893 

0. 1G347545E+0C 

BDl 

2,  4,0.1098) 

2.970 

3.725 

9.893 

0 . 380 1 4066E+00 

BDl 

2,  5,0.109b) 

2.067 

4.292 

9.893 

0 .57774431 E+00 

BDl 

2,  6,0.1098) 

1  .060 

4.645 

9.893 

0. 74637699E  +  00 

BDl 

2,  7,0.109b) 

-0.01,0 

4 . 764 

9.893 

0 . 87758380E  +00 

BDi 

2,  8,0.1098) 

-1  .060 

4.645 

9.893 

0 . 96478379E+00 

BDl 

2,  9,0.1098) 

-2.067 

4 . 292 

9.893 

0 . 1  0036049E+01 

BDl 

2,10,0. 1098) 

-2.970 

3.725 

9.893 

0 . 992 1 0227E+00 

BDI 

2.11,0.109b) 

-3.725 

2.970 

9.893 

0. 93085 104E +00 

BDl 

2,12,0.1 098 ) 

-4.292 

2.067 

9.893 

0.8229221 1 E+00 

BDI 

2,13,0.1 09b ) 

-4.645 

1 .060 

9.893 

0 . 67372888E  +  00 

BDI 

2,14,0.1098) 

-4 . 764 

-0.000 

9.893 

0 . 49075079E+00 

(AMP-M) 

( li  ER  OKU) 

0.000 
0 .  (Jdii 
1  .000 


IESLA 
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BfM  t  t , 12.0. 1098) 
BDf 1 1 , 13,0. 1098) 
BDM  1  , 14,0. 1098) 
B0U1  ,  IS, 0.1038) 
80(11, 16,0.1098) 
BO <11 ,17.0.1098) 
BO ( 11 ,18,0.1098) 
BIH  1  1 ,13.0.1098) 
BD ( 1 1 ,20,0.1 038 ) 
BO ( 1 1 ,21.0.1098) 
BD< 1 1 ,22.0. 1098) 
80(11,23,0.1098) 
80(11,20.0.1098) 
80(11 ,29,0.1 098 ) 
80(11,26.0. 1098) 
BO ( 1 1 .27,0.1098) 
BO (12,  0,0.1098) 
80(12.  1.0.1098) 
BD<  12.  2,0.1098) 
80(12.  3.0.1098) 
DO (12,  4,0.1098) 
BDM2.  5,0.1098) 

80(12,  8.0. 1098) 
80(12.  7.0.1098) 
80(12,  8,0.1098) 
B0M2.  9.0.1096) 
BEX  12.10.0.1098) 
80(12.11  .0.1096) 
BD  (  12, 12.0. 1098) 
80(12.13,0.1 096 1 
80(12.14.0.1098) 
80(12.15.0.1096) 
80(12.16.0.1098) 
BD< 12. 17. 0.1096) 
BD( 12. 18.0. 1098) 
BO  M2. 19.0. 1096) 
BD( 12.20.0. 1098) 
80(12.21 .0.1096) 
80(12.22.0.1098) 
BD ( 1 2 . 23 .0.1 098 ) 
BD ( 12,24.0. 1098) 
80(12.25.0. 1096) 
BO ( 12.26.0. 1098) 
BD< 12.27.0. 1096) 


-6. 1 68 

2.970 

-8.535 

-6.674 

1  . 623 

-6.536 

-6.846 

-0.000 

-3.536 

-6.674 

-1.523 

-6.536 

-6. 168 

-2.970 

-8.536 

-5.352 

-4.266 

-6 . 536 

-4.268 

-6.352 

-8 . 606 

-2.976 

~6.1 06 

-6.536 

-1 .523 

-6.6/4 

-8 . St>6 

C .  006 

-6.846 

-  ,.536 

1  .523 

-6.674 

-8.535 

2.970 

-6.166 

-6 . 586 

4.268 

-6.352 

-8.536 

5.352 

-4.268 

-8 . 536 

6 .  1 68 

-2.970 

-8.586 

6.674 

-1 .623 

-8.536 

4.764 

0.000 

-9.393 

4.645 

1 .066 

-0.893 

4.292 

2.067 

-9.893 

2.726 

2.976 

-9.395 

2.970 

3.726 

-9.893 

2.067 

4 . 292 

-9.395 

1  .060 

4.645 

-9.393 

-0.060 

4  .  764 

-9.893 

-  1  . 060 

4. 645 

-3.393 

-2.06 7 

4 . 292 

-  9 . 89  ? 

-2.970 

3.725 

-9.893 

-3  725  2.970  -9.303 

-4.292  2.067  -9.393 

-4.645  1.060  -9.653 

-4.764  -0.000  -9.893 

-4.645  -1.060  -9.303 

-4.292  -2.067  -9.393 

-3.725  -2.970  -9.603 

-2.970  -3.725  -9.893 

-2.067  -4.202  -9.353 

-1.060  -4.645  -9.893 

O.OCO  -4.764  -9.303 

1.060  -4.645  -9.893 

2.067  -4.292  -9.803 

2.970  -3,725  -3.693 

3.725  -2.970  -9.35? 

4.292  -2.067  -9.393 

4.645  -1.060  -5.30? 


80(13.  0. (1.10981 
BD (13.  1.0.1098) 
80(1 3'.  2,0.1098) 
BD ( 13.  5.0.1098) 
BD (  13.  4.0.1088) 
BD ( 13.  5.0.1098) 
BO ( 13.  6.0. 1098) 
80(13.  7,0.1098) 
80(13.  8.0.1098) 
B 0 ( 13.  9,0.1093) 
BD( 13, 10,0. 1098) 
BO ( 13.1 1 ,0. 1098) 
00(13.12,0. 1098) 
BO  (1 3 . 1 5 . 0 . 1 093 ) 
BD( 13. 14.0. 1098) 
80(13.15.0. 1093) 
80(13.16,0. 1093) 
BO ( 13. 1 7.0.1093) 
80(13,18,0.1093) 
BD ( 13,19.0.1093) 
BD< 13,20,0. 1098) 
80(13.21 ,0.1098) 


2.443 
2.332 
2.201 
1 .910 
I  .  523 
1  . 060 
0.544 
-0.000 
-0.544 
-1 .060 

-i im 
-2.201 
-2.382 
-2.443 
-2.382 
-2.201 
-1 .910 
-1 .523 
-1 .060 
-0.544 
C.OOO 


0 

000 

-in. 

706 

6 

544 

-16. 

706 

1 

060 

-10. 

7  06 

1 

525 

-16. 

706 

1 

910 

-10. 

7  06 

201 

-16. 

7  06 

'■> 

OOO 

1 

-10. 

705 

-• 

,445 

-  1  0  . 

706 

'1 

,382 

-10. 

706 

.201 

-16. 

7  06 

1 

.910 

-  in. 

706 

1 

.523 

-16, 

706 

1 

.  060 

-id. 

706 

C 

.544 

-16 

706 

-0 

■  QUO 

-10 

706 

-0 

.544 

-16 

706 

-1 

.060 

-  10 

705 

-1 

.525 

-16 

706 

-1 

.910 

-  10 

706 

-2 

.201 

-16 

706 

-2 

.  382 

-  HI 

706 

”  cl 

.44  1, 

-16 

.  706 

0.  14  7801608-01 
0.721  86  76/1.-02 
-0.2930!  385C- Hi. 
-0. 721 924336-08 
-0. 147307716-01 
-0.2295r./(l9f--01. 

-0 . 31 90  7S99L-0 1 
-0.4 16661306 -01 
-0.520763926-01 
-0.626364356-01 
-0. 72291 7326-01 
-0.7919571 36  -  0 1 
-0. 807136306-01 
-0 . 7397645/6 -0 1 
-0.572436496  0  1 
-0.513443506 -01 

0 .  on on oo 006  +  00 
0.22 7 946 72’t-  O’ 

0 . 422673 7  46  -  0  1 
0 .562953576-01 
0 . 639.;  1  7046-01 
0.6561 1 1246  -01 

0  .  r.2’6.'  1  89  1  r.  -0  1 
0 .5651,89 1  16  -  01 
0 . 4862  1  9  1  '8c  -  0  1 
0.4005 19x66  -01 
0.313637356- 01 
0 . 229603216-0 1 
0 . 1 5006 ( 586  -  0 1 
6. 7 39/21.. '16  •0_ 
-0,2395 5 3236 -06 
-0  .  753/651  '6  -0.. 
-0 . 160086065-0! 
-0 . 2?9(,  1  1  166-0  1 
-0. 31 364277c -0 ! 
-0 . 400525366-01 
-  0 . 4863258 1 E -  0 1 
-■(1.5(556444/6-01 
-0.626223666-0! 

-n.6561  124  1.6. -01 

-0.639314066-01 
-0.562973206-01 
-0 . 4226.685 36 -0  1 
-0.227  , 0,6  -o' 

0 . 0001)110006+00 

0 . 15742  14  >1  -1)1 

0.239406666-01 

0.41  ?6  .6(1(1!  -01 
0 . 492280.11.6  0  i 
0.6331 :  .3391  0  1 
0 , 6  Jo  7564 41  -  0  1 
0.51  4 4‘ '. ' 9 (  -i'1 
0 . 41-.  7284  1  71  -'ll 
0 .40  .(■."•"/ 46  cl 
0  .  j'70  I  (2.746  u  I 
0. 25051  '3296  -'ll 
0.  1  6797.4i.3l  -  0  1 
0 . 1.4  1  4'  i.'lbSi  -ll7 
-fl.  2  7404  I!  66  oi, 
-0 .54165115  It  07 
-0.167940181  01 

-0 .26115.  00  21  -di 

-0 . 780  I  70  7)1  0  1 
-0 . 411  i3536  11 1 

-0 . 46 7286.266  -n  1 
-0 .51  461.5066  'll 


It 

B 

u 

2  *  * 

n 


BD  ( 

9.  2,0.1098) 

8.313 

4.292 

-4.764 

BD  ( 

9,  3,0.1098) 

7.734 

6.168 

-4.764 

BD  ( 

9,  4.0.1098) 

6.163 

7.734 

-4.764 

BD  ( 

9.  5,0.1098) 

4 . 292 

8.913 

-4.764 

BD  ( 

9.  6.0.1098) 

2.201 

9.645 

-4.764 

BD  ( 

9,  7,0.1098) 

-c.ooc 

9.395 

-4.764 

BD  C 

9.  8,0.1098) 

-2.20! 

9.645 

-4.764 

BD  ( 

9,  9.0.1098) 

-4.292 

8.913 

-4.764 

BD  ( 

9.10,0.1098) 

-6.168 

7.734 

-4.764 

BD  ( 

9,11.0.1093) 

-7 .  734 

6 . 1 66 

-4.764 

BD  ( 

9.12.0.1098) 

-8.913 

4.292 

-4.764 

BD  ( 

9.13.0.1098) 

-9.645 

2.201 

-4.764 

BD  ( 

9,14.0. 1098) 

-3.893 

-0.0  '0 

-4.764 

BD  ( 

9.15,0.1 038 ) 

-3.645 

-2.2  )1 

-4.764 

BD( 

9.16,0. 1098) 

-3.913 

-4.292 

-  4 . 764 

BD  C 

9.17,0.1 098 ) 

-7.734 

-6.168 

-4.764 

BD( 

9,18,0.1 098 ) 

-6. 168 

-7.7 34 

-4 . 764 

BD  ( 

9.19,0.1 098 ) 

-4.292 

—  o .  3  i  j 

-4. 764 

BD  ( 

9.20,0. 1038) 

-2.20  1 

-9.645 

-4 . 764 

BD  I 

9.21,0.1098) 

C.OOC 

-3.833 

-4.764 

BDl 

3.22,0. 1098) 

2.201 

-9.645 

-4.764 

BD  f 

9.23,0.1 033) 

4.292 

-6.315 

-4.  7 64 

BDl 

9.24,0. 1098) 

6.163 

-7.734 

-4 . 764 

BDl 

9,25.0. 1038) 

7.764 

-6.163 

-4.764 

BDl 

9.26,0.1098) 

8.313 

-4.292 

-4.764 

BDl 

9.27,0.1098) 

9.645 

-2.2C1 

-4.764 

BDdO,  0  .  j  .  1038) 

8.534 

0.000 

-6.846 

BDdO,  1,0.1098) 

8.363 

1.910 

-6.846 

BDdO,  2,0.1038) 

7.734 

3.725 

-6.846 

BD c 10.  3.0.109S) 
BD ( 10,  4,0.1098) 
BD <  10.  5,0.1098) 
BD ( 1 0 ,  6,0.1098) 
BD ( 1 0 ,  7,0.1098) 
BDdO,  8,0.1098) 
BD ( 1 0 ,  9.0.1098) 
BD  <  10, 10,0. 1098) 
BDdO.1  1  ,0.1098) 
BD(10, 12,0.1098) 
BD (10. 13,0. 1098) 
BDdO,  14,0. 1098) 
BD(10, 15.0.1098) 
BD ( 10, 16,0. 1098) 
BDdO, 1  7,0. 1098) 


6.712 

5.352 

3.725 


5 . 352 
6 !  7 1*2 
7.734 


■  6 . 846 
-6.846 
-6.846 


BD (  10, 18.0 
BDdO.  19,0 
BD( 10.20.0 
BDdO, 21 .0 
BDdO,  22,0 
BDdO, 23,0 
BDdO, 24,0 
BDdO,  25.0 
BDdO, 26,0 
BDdO 
BD  (  1  1 
BDd  1 


BD  (  1  1 
BDd  1 
BD  (  1  1 
BDd  1 
BDd  1 
BD  d  1 
BDd  1 
BDd  1 
BDd  1  . 

BDd  1  . 1 1  .0. 


?7.0. 
0,0. 
1  .0. 
2,0. 
3,0. 
4.0. 
5.0. 
6,0. 
7,0. 
8.0. 
9.0. 
10.0 


.1098) 
.109.3) 
.1098) 
.1  093) 
.1098) 
.1  093) 
.1098) 
.1098) 
.1098) 
.1098) 
.1098) 
.1038) 
.1098) 
.1098) 
.1098) 
1038) 
1098) 
1  098) 
1098) 
1098) 
1098) 
1098) 


1.910 

3.369 

-6.846 

-o.oco 

3.584 

-6.846 

-1.910 

8.369 

-6.846 

-3.725 

7.7  34 

-6.846 

-5.352 

6.712 

-6.846 

-6.712 

5.352 

-6.846 

-7.734 

3.725 

-6.846 

-8.369 

1.910 

-6.846 

-8.534 

-0.000 

-6.846 

-3.369 

-1 .910 

-6.846 

-7.734 

-3.725 

-6.846 

-6.712 

-5.352 

-6.346 

-5.352 

-6.712 

-6.846 

-3.725 

-7.734 

-6.846 

-1.310 

-8.369 

-6.846 

C.OOC 

-8.584 

-6.846 

1  .910 

-8.369 

-6.846 

3.725 

-7.734 

-6.846 

5.352 

-6.712 

-6.846 

6.712 

-5.352 

-6.846 

7.734 

-3.725 

-6.846 

8.369 

-1 .910 

-6.846 

6.846 

0.000 

-8.535 

6.674 

1  . 523 

-6.535 

6.168 

2.970 

-8.535 

5.352 

4 . 26,9 

-8.585 

4.268 

5.352 

-8.535 

2.970 

6 . 1  68 

-8.535 

1  _  C  TO 

6.674 

-3.585 

-c.ooc 

6.846 

-8.635 

-1 .523 

6.674 

-3.535 

-2.97C 

6.163 

-8.585 

-  4 . 268 

5 . 352 

-3.535 

-5.352 

4.268 

-8 . 585 

-0.32168325E-01 
-0.24101 S44E-0 1 
-0. 14826197E-01 
-0.48073940E-02 
0.54524504E-02 
0. 1S438925E-01 
0 . 24651 207E-0 1 
0.32627348E-01 
0.38967434E-01 
0.43353509E-01 
0.45565665E-01 
0.45492953E-01 
'  0.431 39022E-0 1 
0 . 3362 1391 E -01 
0.321 63072E-0 1 
0 . 24 1 0 1 272E-0 1 
n. I 4825344E-0 1 
0.48071072E-02 
-0.54527447E-02 
-0. 1543921 6E-01 
-0.24651498E-01 
-0.3262/620E-01 

-0.38967643E-0I 
-0.4335361 7E-01 
-0 . 4556569 1 E-0 1 
-0.45492906E-01 
-0.29663622E-01 
-0 .2661 91 36E-01 
-0.22239853E-01 
-0. 16745351E-01 
-0.104 1 ! 192E-01 
-0.35549663E-02 
0 . 347952 1 5E-02 
0 . 1 0329562E-01 
0. 1S681086E-01 
0.221851 75E-U1 
0.26578754E-01 
0.29638540E-01 
0.31212129E-01 
0 . 31 220622E-01 
0 . 29663540E-0 1 
0 . 2661 9025E-0 1 
0.22239670E-01 
0. 16745169E-01 

0.  10410976E-01 
0.35547570E-02 
-0 . 34797254E -02 
-0. 10339763E-01 
-0. 16681295E-01 
-0.22136365E-01 
-0 . 2657888 1 E-0 1 
-0.2963S622E-01 
-0.312121 55E- 0 1 
-0.31 22C535E- 0 1 
-0. 1 983S476E-0 1 
-0.1 7822828E-01 
-0. 1 49 1 6483E-0 1 
-0.1 1 262 1 64E-0 1 
-0. 70431 009E-02 
-0 . 247067 42E -  02 
0 . 2225256  1 E  -02 
0.6809S232E-02 
0. 1 1052388E-0! 

0 . 1 4741 715E-01 
0  .  1  7C.31  3ME-01 

0. 19/538406-01 
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I 


fc 


BD< 11,12.0. 1098) 

-6.168 

2.970 

_  D  COE 

'J  .  JU  J 

BD( 11 .13.0.1098) 

-6.674 

1  .523 

-S.565 

BD( 1 1 . 14.0. 1098) 

-6.846 

-0.000 

-8.585 

ED ( 1 1 .15.0.1098) 

-6.6/4 

- 1  . 523 

-8.565 

BD( 1 1 . 16.0. 1098) 

-6 . 168 

-2.970 

-8.585 

BD ( 1 1  .17,0.1098) 

-5.352 

-  4 . 268 

-8.565 

BD( 1 1 . 13,0. 1098) 

-4.268 

-5.352 

-8.585 

BD ( 1 1 .19,0.1098) 

-2.970 

-6.168 

-8.565 

BD( 1 1,20.0.1098) 

-1.523 

-6.574 

-8.585 

BP( 1  1 ,21  .0.1098) 

O.OCO 

-6.846 

-8.565 

BD  f  1  1,22,0. 1098) 

1  .523 

-6.674 

-8.585 

BDC11 ,23.0.1098) 

2.970 

-6.168 

-8.585 

BD( 1 1 ,24,0. 1098) 

4.268 

-5.352 

-8.585 

BD (11 ,25,0.1098) 

5.352 

-4.268 

-8.565 

BD ( 1 1 ,26,0. 1098) 

6.168 

-2.970 

-8.585 

BD( 11 .27.0.1098) 

6.674 

-1 .523 

-8.565 

BD ( 12,  0.0.1098) 

4.764 

0.000 

-9.893 

BD ( 1 2 ,  1,0.1098) 

4.645 

1  .060 

-9.893 

BD ( 12,  2,0.1098) 

4.292 

2.067 

-9.893 

BD ( 1 2 .  3,0.1098) 

3.725 

2.970 

-9 . 893 

BD ( 1 2 ,  4,0.1098) 

2.970 

3.725 

-'9.893 

BD ( 1 2 ,  5,0.1098) 

2.067 

4.292 

-9.893 

BD ( 12,  6.0.1098) 

1.060 

4.645 

-9.893 

BD ( 1 2 ,  7,0.1098) 

-c.ooc 

4.764 

-9.893 

BD ( 12,  8,0.1098) 

-1.060 

4.645 

-9.893 

BD ( 1 2 ,  9,0.1098) 

-2.067 

4 . 292 

-9.893 

BD( 12, 10,0. 1098) 

-2.970 

3.725 

-9.893 

BD ( 1 2 . 1  1  ,0.1098) 

-3.725 

2.97C 

-9.893 

BD( 12. 12,0.1098) 

-4.292 

2.067 

-9.893 

BD ( 1 2 , 13,0.1098) 

-4.645 

1 .060 

-9.893 

BD( 12, 14,0. 1098) 

-4.764 

-0.000 

-9.893 

BD ( 1 2 . 15, 0.1098) 

-4.645 

-1 .06C 

-S.893 

BD( 12, 16,0. 1098) 

-4.292 

-2.067 

-9.S93 

BD ( 1 2 . 17, 0.1098) 

-3.725 

-2.97C 

-9.893 

BD  C 12, 18.0.1098) 

-2.970 

-3.725 

-9.893 

BD ( 1 2 . 19, 0.1098) 

-2.067 

-4.232 

-9.393 

BD  C 12.20.0. 1098) 

-1.060 

-4.645 

-9.893 

BD  c 12.21  ,0.1098) 

C.OOC 

-4.764 

-9.893 

BD ( 12.22.0. 1098) 

1  .060 

-4.645 

-9.893 

BD (12. 23, 0.1093) 

2.067 

-4.292 

-9.893 

BD ( 12.24,0.1098) 

2.970 

-3.725 

-9.893 

BD(12. 25, 0.1098) 

3.725 

-2.970 

-9.893 

BD  < 12.26,0.1098) 

4.292 

-2.067 

-9.393 

BD< 12,27,0.1093) 

4.645 

-1 . 06C 

-9.396 

BD ( 1 3 ,  0,0.1098) 

2.443 

0.000 

-10.705 

BD ( 1 3 ,  1,0.1098) 

2.382 

0.544 

-10.7C5 

BD ( 13.  2.0.1098) 

2.201 

1  .060 

-10.705 

BD ( 13.  3.0. 1098) 

1  .910 

1  .523 

-10.7C5 

BD ( 13,  4,0.1098) 

1  .523 

1.910 

-10.705 

BD ( 1 3 .  5,0.1098) 

1.060 

2.2C1 

-10.7C5 

BD ( 13,  6.0.1098) 

0.544 

2.382 

-10.705 

BD ( 1 3 ,  7,0.1098) 

-O.OCO 

2.443 

-10.7C5 

BD ( 1 3 ,  8.0.1098) 

-0.544 

2.382 

-10.705 

BD ( 1 3 ,  9,0.1098) 

-1.060 

2.2C1 

-10.7C5 

BD ( 13. 10,0. 1098) 

-1.523 

1.910 

-10.705 

BD ( 13. 1 1  ,0.1098) 

-1 .910 

1  .523 

-10.7C5 

BD( 13. 12,0. 1098) 

-2.201 

1.060 

-10.705 

BD ( 13. 13.0. 1096) 

-2.362 

0.544 

-10. 7C5 

BD< 13.14,0.1098) 

-2.443 

-0.080 

-10.705 

BD ( 1 3 , 15. 0.1098) 

-2.382 

-0.544 

-10.7C5 

BD  < 13, 16,0. 1098) 

-2.201 

-1.060 

-10.705 

BD ( 13. 17, 0.1096) 

-1 .910 

-1  .523 

-10.  7C5 

BD( 13. 18,0. 1098) 

-1.523 

-1.910 

-10.705 

BD ( 1 3 , 19, 0.1096) 

-1.060 

-2.2C1 

-10.7C5 

BD  C 13,20,0. 1098) 

-0.544 

-2.382 

-10.705 

BD ( 13,21 ,0.1098) 

O.OCO 

-2.443 

-10. 7  C  5 

0.20825738E-01 
0 . 203S3449E-01 
0. T  9S3S^3 IE-01 
0.17822765E-01 
0. 1 49 1 6368E-0 1 
0.11 2B204AE-01 
0 . 70429705E-02 
0.247C73SSE-02 
-0.222S392SE-02 
-0.68099536E-02 
-0.11 0B3022E-0 1 
-0.14741 842E-01 
-0.1 7691 407E-0 1 
-0.1975388SE-01 
-0.2082S807E-01 
-0.208S3437E-01 
-0. 12203962E-01 
-0.10973342E-01 
-0.91 924742E-02 
-0.69SCSS04E-02' 
-  0 . 436O29S8E-02 
-0.1 551 2977E-02 

0. 13354912E-02 
0 . 41 SS3229E-02 
0 . 676S7738E-02 
0 . 903891 40E-02 
0. 108S7802E-01 
0.121 32239E-01 
0. 12798298E-01 
0. 12322G09E-01 
0. 1220393GE-01 
0 . 1 0973297E-01 
0.91924034E-02 
0.6950579GE-02 
0.43G02139E-02 
0. 1S512137E-02 
-0. 133558G0E-02 
-0.4155401  IE-02 
-0.67GG8483E-02 
-0 . 90389922E-02 
-0. 10857854E-01 
-0.121 322S5E-01 
-0.127983  9E-01 
-0. 12822602E-01 

-0.581 93840E-02 
-0 . 5234334GE-02 
-0 . 438681 99E-02 
-0.331 9330SE-02 
-0 . 20853961 E-02 
-0 . 74688881 E-03 
0 . S2907068E-03 
0 . 1 9734891 E-02 
0.321 89409E-02 
0.4302986CE-02 
0.51 712580E-02 
0.578C2238E-02 
0.60993396E-02 
0.611 26091 E-02 
0.58! 9369 1 E -  02 
0 . 523431 23E-02 
0.43867864E-02 
0.331 92968E-02 
0 . 20853553E-02 
0 . 74664294E-03 
-0.6291 1375E-03 
-0. 19735277E-02 
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BD( 13.22.0. 1098) 
BD ( 1 3 , 23, 0.1038) 
BD ( 1 3 , 24 , 0 . 1 098 ) 
8DM3, 25,0.  1093) 
BD ( 13,26,0. 1098) 
80(13,27,0.1098) 
BD ( 1 4 ,  0.0.1098) 
END  OF  DATA 


0 

.544 

-2, 

.382 

-10. 

.705 

1 

.  06C 

-2 , 

.201 

-1C. 

,705 

1 

.523 

-1  . 

.910 

-10. 

,705 

1 

.91  C 

-1  . 

.523 

-10. 

,705 

2 

.201 

-1  , 

,060 

-10. 

,705 

2 

.382 

-C 

.544 

-1C, 

.  705 

-0 

.000 

0 

.000 

-10, 

.980 

-0 . 32 1 89772E -  02 
-0 . 4303C 1 95E-02 
-0.5171 2804E-02 
-0 . 57802387E-02 
-0.60993433E-02 
-0.G1126053E-02 
-0.831471 78E-08 
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Data  for 


DIF'3LE*1  DIP0LE*2  DIP0LE*3 

moment  A.11E-U4  0.1  IE-04  O.OOE+OO  (AMP*M) 

LOCATM  A .  0  7  A  0.070  0.000  ( li  FR  ORN) 

trans  0.711,-,  0. aao-h.opo  0.000-0. 906-0.423  1.000  0.000  0.000 

MATRIX  U.OOO  1.000  U.000  1.000  O.OGO  0.000  0.000  1.000  0.000 

0.023  O.OIIO  0.9AF,  0.000-0.023  0.906  0.000  0.000  1.000 

THE  HEAS  SPHERE  RADIUS  IS  0.1098  M 

SPHER  DATA  POIIIIS:  1ULTA-  I  »a  .2244;  PHI  *F -‘0.2244 

BD(T.F.R)  X  Y  Z  FLUX  DENS,  XEE-11  TESLA 


[!D( 

0,  0,0.1  trio ) 

0 . 000 

0.000 

10.980  . 

.  --0. 12354885E-05 

BD( 

1,  0,0.1098) 

2 . 44  J 

G.OOG 

1G. 705 

0.781 14372E+00 

BD( 

1,  1,0.1098) 

2 . 382 

0.544 

10.705 

0.11 764374E+0 1 

BO  ( 

1  ,  2,0.1098) 

2.201 

1  .060 

1G.705 

0 . 1 4506912E+01 

BD  f 

1,  3, U. 1098) 

1.910 

1  . 523 

10.705 

0. 1S39555SE+01 

B!>  < 

1,  4.  A. 1098) 

1 . 523 

1  .910 

1G. 705 

0. 14S97769E+01 

BD  ( 

1.  S. 0.1098) 

1  .060 

2.201 

10.705 

0. 12703056E+01 

BD  ( 

1  .  A  ,0. 1098) 

0.544 

2.382 

10.705 

0. 10306473E+01 

BD  ( 

1.  /,  0.1 090) 

-G.OOG 

2.443 

1C. 705 

0 . 781 1 3908E+00 

BD  ( 

1,  8, U. 1098) 

-0 . 544 

2.382 

10.705 

0.541 58998E+00 

BD  ( 

1  ,  9,0.1  098) 

-1 .060 

2.201 

10.705 

0.31 708521 E  +  OO 

BD  ( 

1.10,0. 1098) 

-1.523 

1.910 

10.705 

0. 10426S1SE+00 

BD  ( 

1.11  ,0.1099) 

-1 .910 

1  .523 

10.705 

-0.1 04271 89E+00 

BD  ( 

1,12,0.1098) 

-2.201 

1 .060 

10.705 

-0.3171001 7E+00 

BD  ( 

1,13,0.1 098 ) 

-2.382 

G.544 

10.705 

-0.541601 60E+00 

GD< 

1,14,0.1098) 

-2.443 

-0.000 

10.705 

-0.781 1SS35E+00 

BD  ( 

1  .IS,  0.1090) 

-2.382 

-0.544 

10.705 

-0.10306S49E+01 

BD  < 

1,16,0.1098) 

-2.201 

-1.060 

10.705 

-0. 127031 14E+01 

BD  ( 

1,17,0.1098) 

-1 .910 

-1 .523 

10.705 

-0.1 4597826E+0 1 

BD( 

1,18,0.1098) 

-1.523 

-1.910 

10.705 

-0. 15395603E+01 

BD  ( 

1.19,0.1098) 

-1 .060 

-2.201 

1C. 705 

-0.1 4506731 E+0 1 

BD( 

1,20,0. 1098) 

-0.544 

-2.382 

10.705 

-0. 1 1764164E+01 

BD  ( 

1,21,0. 1098) 

O.OOG 

-2.443 

1C.  705 

-0.781 12715E+00 

BD  ( 

1,22,0.1098) 

0.544 

-2.382 

10.705 

-0 . 39578557E  +  00 

BD  ( 

1 ,23,0.1098) 

1  .066 

-2.201 

10.705 

-0.1  3669676E  +  00 

BD( 

1.24,0.1098) 

1  .523 

-1.910 

10.705 

-0.24487235E-01 

BD  ( 

1 ,29,0.1098) 

1  .910 

-1  .523 

10.705 

0 . 24488252E-0 1 

BD  ( 

1,26,0.1098) 

2.20! 

-1.060 

10.705 

0.  13672072E  +  00 

BD  ( 

1  ,27,0.1098) 

2 . 382 

-G.544 

10.705 

0.39582533E+00 

BD  ( 

2.  0.U.1098) 

4.764 

0.000 

9.893 

0 . 822 1 5434E+00 

BD( 

2.  1,0.1096) 

4.645 

1  .060 

9.893 

0 . 21 3321 S9E  +  0 1 

BD  C 

2,  2,0.1098) 

4.292 

2.067 

9.893 

0.27734528E+01 

BD  ( 

2.  3,0.1096) 

3.725 

2.970 

9.893 

0 . 26963927E+0 1 

BD  ( 

2.  4.0.1098) 

2.970 

3.725 

9.893 

0 . 22320309E  +  0  1 

BD  ( 

2,  S, 0.1 096) 

2.067 

4 . 292 

9.893 

0. 168S57SSE+01 

BD( 

2.  6,0.1098) 

1  .060 

4.645 

9.893 

0. 120401 10E+01 

BD  ( 

2.  7.0.1096) 

-O.OGO 
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Appendix  D 


This  section  contains  the  data  collected  from  each  of 
the  five  multipole  conducting  sphere  experiments. 


Experiment  Data  for  Model  H 
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